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TECHNIQUES  FOR  THE  ANALYSIS  OF  EXPENDABLE  BATHYTHERMOGRAPH  (XBT) 
DATA  WITH  APPLICATIONS  TO  THE  SOUTH-WESTERN  PACIFIC  OCEAN 


P.G.  Marshallsay  and  S.M.  Ball 


SUMMARY 

The  expendable  bathythermograph  is  a  convenient  and 
inexpensive  tool  for  the  measurement  of  ocean  temperature 
structure.  It  does  not  measure  salinity,  but  in  many 
areas  of  the  ocean  it  is  reasonable  to  use  standard 
temperature-salinity  (T-S)  curves  to  infer  salinity  from 
temperature  measurements.  T-S  curves  constructed  from 
historical  Nansen  cast  data  are  presented  for  10°  squares 
in  the  South-Western  Pacific  Ocean,  covering  the  area  from 
15°S  to  45 °S,  and  from  150°E  tv.  180°E.  With  a  knowledge 
of  both  temperature  and  salinity  structures  a  wide  range 
of  other  parameters  may  be  estimated,  the  most  important 
of  these  being  density,  dynamic  height  and  sound  speed. 
Relationships  for  the  computation  of  these  and  other 
properties  of  sea  water  are  discussed,  and  computer 
programmes  presented  to  calculate  profiles  of  the  relevant 
parameters  from  XBT  data . 


POSTAL  ADDRESS:  Chief  Superintendent,  Weapons  Systems  Research  Laboratory, 
Box  2151,  GPO,  Adelaide,  South  Australia,  5001. 
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1.  INTRODUCTION 

The  oceanographer  is  frequently  confronted  with  the  need  to  make  extensive 
measurements  of  large-scale  features  in  a  scenario  which  may  undergo 
significant  changes  in  a  relatively  short  span  of  time.  In  this  situation 
the  need  for  near-synoptic  information  covering  the  area  of  interest  is 
obvious,  and  the  delays  inherent  in  the  traditional  techniques  of  making 
Nansen  casts,  or  more  recently  conductivity-temperature-depth  (CTD)  surveys, 
which  involve  prolonged  stoppages  on-station  may  be  highly  undesirable,  if  not 
unacceptable.  Against  this  background  the  expendable  bathythermograph  (XBT) 
has  established  itself  over  the  past  decade  or  so  as  an  invaluable  tool  in  the 
conduct  of  rapid  oceanographic  surveys,  either  supplementing  or  replacing  the 
previously-mentioned  techniques .  The  XBT  briefly  comprises  a  low-cost, 
weighted,  finned  projectile  with  a  thermistor  bead  mounted  in  the  nose, 
together  with  associated  launching  equipment  and  recording  instrumentation. 
The  thermistor  forms  one  arm  of  a  bridge,  into  which  it  is  connected  by 
several  hundred  meters  of  lightweight,  insulated  copper  wire,  which  is  spooled 
out  as  the  projectile  falls.  The  hydrodynamic  characteristics  of  the 
projectile  produce  a  known  rate  of  descent,  so  that  temperature  may  be 
recorded  as  a  continuous  function  of  depth.  The  value  of  the  technique  lies 
in  the  fact  that  measurements  may  be  made  in  transit,  and  a  given  area  may 
thus  be  surveyed  in  a  considerably  lesser  time  than  would  otherwise  be 
possible. 

A  survey  conducted  using  such  techniques  produces  a  near-synoptic  picture  of 
the  temperature  structure  of  the  area  under  investigation.  The  rapidity  and 
relatively  low  cost  with  which  such  a  survey  may  be  performed  makes  the 
optimal  use  of  the  accruing  data  as  a  means  of  inferring  details  of  the 
distribution  of  dynamic  height  and  sound  speed,  together  with  other  properties 
of  the  ocean  waters,  highly  desirable.  Estimation  of  such  properties 
requires  however,  that  the  corresponding  salinity  structures  of  the  areas 
under  study  also  be  known.  It  has  long  been  recognised  that  for  most  of  the 
World's  ocean  waters  a  close  relationship  exists  between  temperature  and 
salinity,  and  a  number  of  investigators  have  examined  the  possibility  of  using 
this  relationship  to  infer  salinity  from  temperature  measurements  made  using 
bathythermographs  (BTs)  and  XBTs(refs.l  to  3).  The  results  of  such  studies 
have  been  sufficiently  en  'ouraging  to  induce  a  large  number  of  workers  to 
adopt  the  use  of  temperature-salinity  (T-S)  curves  derived  either  from  Nansen 
or  CTD  surveys  made  concurrently  with  an  XBT  survey,  or  in  the  absence  of  such 
measurements,  from  historical  data,  as  a  standard  means  of  processing  XBT  data 
to  derive  estimates  of  temperature  and  salinity-dependent  oceanographic 
variables  (see  for  instance,  reference  4).  Such  temperature-salinity 
relationships  tend  to  be  applicable  over  relatively  wide  geographical  areas, 
and,  with  the  exception  of  the  topmost  tens  of  metres  of  the  ocean,  which  are 
subject  to  seasonal  and  shorter  term  influences,  to  be  substantially 
independent  of  depth.  In  a  survey  of  available  historical  data  for  the 
Pacific  Ocean  between  20°S  and  40°N,  considered  as  being  divided  into  10° 
squares,  Emery  and  Wert(ref.5)  found  the  error  in  estimating  dynamic  height 
from  temperature  measurements  augmented  by  salinities  inferred  using  standard 
T-S  relationships  to  exceed  uncertainties  in  the  measurements  accruing  from 
other  causes  in  only  about  25%  of  the  areas  considered.  These  they 
identified  as  being  transition  regions  between  different  water  masses. 
Flierl(ref .6)  further  considered  the  problem  of  reducing  XBT  measurements  made 
in  transition  zones  such  as  may  be  found  on  the  boundaries  of  eddies,  and 
derived  a  procedure  for  the  processing  of  such  data  in  conjunction  with 
supplementary  topical  CTD  measurements. 
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In  the  present  study  we  address  ourselves  to  the  case  where  the  use  of 
standard  T-S  curves  may  be  regarded  as  providing  the  desired  accuracy  and 
reliability,  and  discuss  the  derivation  of  profiles  of  properties  of  interest 
using  such  methods.  A  series  of  standard  T-S  curves  for  the  South-West 
Pacific  Ocean  between  15°S  and  45°S,  and  between  150°E  and  180°E  are  developed 
from  historical  data,  and  their  implementation  as  a  basis  for  a  computational 
scheme  considered.  Relationships  for  the  computation  of  dynamic  height  and  a 
variety  of  seawater  properties  dependent  upon  temperature,  salinity  and 
pressure  are  then  presented.  Finally,  a  series  of  computational  routines  to 
process  XBT  data  to  produce  profiles  of  the  previously  introduced  properties 
are  listed  and  described. 


2.  T-S  CURVES  FOR  THE  SOUTH-WESTERN  PACIFIC  OCEAN 

For  the  analysis  of  XBT  data  from  the  Tasman  Sea,  the  area  from  15°S  to  45°S, 
and  from  150°E  to  180°E  was  divided  into  10°  squares.  The  location  and 
labelling  of  the  squares  is  shown  in  figure  1.  Representative  standard  T-S 
curves  were  then  constructed  using  historic  Nansen  cast  data  originating  from 
CS.vO  and  RANRL  cruises.  The  distribution  of  the  available  data  over  the 
area  of  interest  varied  enormously;  the  region  immediately  beyond  Sydney 
Harbour  was  represented  by  an  extremely  heavy  concentration  of  data,  while  the 
available  data  for  some  of  the  more  remote  squares  was,  by  comparison,  sparse. 
Approximately  twelve  Nansen  cast  profiles  were  chosen  from  each  square  to  form 
a  data  base  for  further  processing.  The  representative  Nansen  casts  were 
selected  so  as  to  be  more-or-less  evenly  distributed  over  their  respective 
areas.  In  addition,  since  it  was  desired  to  pool  all  selected  data  for  a 
particular  square  without  regard  to  season,  wherever  possible  the  casts  were 
selected  so  that  approximately  three  profiles  from  each  season  were  available 
for  a  given  square.  Details  of  the  Nansen  casts  selected  are  reproduced  in 
Appendix  I,  and  the  locations  of  the  casts  shown  on  figure  1. 

T-S  curves  constructed  using  the  pooled  Nansen  cast  data  points  for  each  area 
are  reproduced,  without  identification  of  the  individual  soundings,  in 
figure  2.  Data  points  from  within  the  mixed  layer  had  previously  been  edited 
out  of  the  data  base  by  hand.  The  comparatively  small  size  of  our  data  base 
precluded  the  possibility  of  performing  any  pre-processing  quality  control  of 
the  type  used  by  Emery  and  Wert (ref. 5),  and  by  Pearce  and  Hamon(ref . 7) . 
Nevertheless,  the  data  points  selected  are  seen  to  define  a  series  of  T-S 
curves  which  are  'tight'  in  the  sense  of  reference  2.  As  expected,  scatter 
is  largest  at  the  highest  temperatures,  corresponding  to  the  upper  layers  of 
the  ocean. 

For  computational  purposes  it  is  desired  that  the  T-S  curves  of  figure  2  be 
represented  by  a  series  of  functions  in  a  manner  which  is  in  some  sense 
optimal.  In  the  first  instance  it  was  attempted  to  fit  fourth-order 
polynomials  to  the  data.  Here,  as  in  reference  7,  it  was  found  impossible  to 
faithfully  represent  the  T-S  curves  for  this  geographical  area  by  low-order 
polynomials.  The  problems  become  most  evident  in  the  case  of  those  areas 
such  as  B1  and  Cl  in  which  the  T-S  curves  exhibit  a  curious  'hook'  close  to 
the  minimum  temperature.  In  attempting  to  provide  an  optimal  fit  through 
this  section  of  the  data,  the  polynomial  approximations  were  found  to 
consistently  overestimate  the  deepwater  salinity  minima.  In  order  to 
overcome  this  problem,  Pearce  and  Hamon  omitted  all  points  representing 
temperatures  lower  than  2.5°C  from  the  polynomial  fitting  process.  We  chose 
instead  to  retain  all  data  points,  and  to  represent  the  curves  using  cubic 
splines  fitted  in  the  least-squares  sense.  This  was  done  using  subroutine 
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VC03A  from  the  Harwell  Subroutine  Library(ref .8) .  This  employs  a  least- 
squares  cubic-spline  fitting  algorithm  which  selects  its  own  level  of 
smoothing,  and  chooses  its  own  knot  points  to  span  the  domain  over  which  the 
function  is  defined  by  the  available  data.  The  fitted  curves  correspond  to 
the  solid  lines  in  figure  2,  and  are  seen  to  comprise  an  excellent 
representation  of  the  data,  particularly  at  the  low  temperature  end  of  the 
domain.  The  cubic  spline  coefficients  are  listed  in  Appendix  II. 

Further  appreciation  for  the  relevance  of  the  T-S  relationships  may  be  gained 
by  considering  the  temperature  structure  of  an  ocean  water  column.  For  this 
purpose  the  Nansen  cast  temperature  data  have  been  plotted  against  depth  for 
each  area  in  figure  3.  These  provide  some  indication  of  the  depth 
appropriate  to  a  specific  point  on  the  T-S  curves.  It  should  be  borne  in 
mind  however  that  what  is  postulated  is  not  an  unvarying  temperature-depth 
relationship,  but  a  unique  temperature-salinity  relationship.  Thus,  although 
the  water  at  a  given  depth  within  a  particular  10°  square  may  show 
considerable  variation,  the  salinity  corresponding  to  a  given  temperature  may, 
to  a  close  approximation,  be  regarded  as  being  specified  by  the  appropriate  T- 
S  curve  without  regard  to  the  depth  or  time  of  the  year  at  which  the  sample 
was  taken.  In  general,  these  comments  are  true  for  all  but  the  uppermost 
layers  of  the  ocean.  For  instance,  the  T-S  relationship  for  area  A1  is  seen 
to  be  very  'tight'  indeed  for  temperatures  below  about  20°C.  Higher 
temperatures  than  this  are,  on  the  whole,  confined  to  the  uppermost  200  m  of 
the  ocean,  although  it  will  be  noted  that  at  say,  500  m,  temperature 
fluctuations  over  a  3°  span  are  evident,  within  a  temperature  range 
corresponding  to  a  'tight'  portion  of  the  T-S  curve. 

The  confidence  with  which  the  standard  T-S  curves  may  be  used  to  reduce  data 
from  the  uppermost  layers  of  the  ocean  will  depend  in  part  upon  the  use  which 
is  to  be  made  of  the  results.  Thus,  A.F.  Pearce  (personal  communication) 
suggests  that,  if  the  results  are  to  be  used  to  estimate  dynamic  height,  then 
there  is  little  point  in  attempting  to  correct  the  results  for  surface 
effects,  since  the  uppermost  layers  of  the  ocean  contribute  little  to  the 
dynamic  height.  Andrews ( ref . 3)  further  regards  the  salinity  obtained  from  a 
T-S  relationship  as  comprising  a  second-order  correction  to  the  estimates  of 
the  physical  properties  of  seawater.  He  does  however,  propose  a  correction 
technique  for  surface  effects  which  may  be  used  if  surface  salinity  has  been 
measured,  and  which  comprises  a  technique  for  interpolation  between  the 
surface  salinity  and  the  deep-sea  salinity  inferred  from  a  standard  T-S  curve. 
Thus,  if  the  depth  D  is  selected  as  the  depth  at  which  the  thermal  gradient 
evaluated  over  a  10  m  interval  of  depth  is  a  maximum,  and  L  is  the  depth  of 
the  mixed  layer,  defined  as  the  depth  at  which  the  temperature  has  fallen  to 
98.5%  of  the  average  value  above  that  point,  then  Andrews  proposes  that  the 
salinity  in  the  upper  layers  of  the  ocean  may  be  found  according  to  the 
equations 


S(z)  =  So, 


O^zSL,  (1(a)) 


S(z)  =  (z  -  L) 


S(T(z) )  -  So 

1.5D  -  L 

m 


+  S0, 


LSzS1.5D  , 


(Kb)) 
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and 


S(z)  =  S(T(z) ) , 


zgl.SD  .  (1(c)) 

m 


where  S0  is  the  surface  salinity,  which  is  assumed  to  remain  constant  through 
the  mixed  layer. 

Using  the  measurements  of  temperature  and  salinity  from  Nansen  cast  data,  the 
corresponding  one-atmosphere  densities  have  also  been  computed  using  the 
equation  of  state  of  seawater  as  derived  by  Chen  and  Millero(ref .9) .  The  data 
are  plotted  by  area  in  figure  4,  where  the  one-atmosphere  density  is  expressed 
in  the  familiar  sigma-t  form( ref . 10 , 1 1 ) : 


Ot  =  1000(p°  -  1)  (g. cm. -3 . 10-3)  (2) 

where  p°  is  the  one-atmosphere  expressed  in  g.cm-3.  Andrews(ref . 12)  has 
proposed  that  the  density  structure  in  the  Tasman  Sea  may  be  approximated  by  a 
relationship  of  the  form 


ctt  =  p0  +  dp[ 1  -  1/(1  -  z/D)] 


(3) 


where  the  parameters  p0,  dp  and  1)  remain  to  be  specified  for  a  given  area. 
Optimum  values  for  these  parameters  were  computed  using  a  nonlinear  least- 
squares  function  fitting  subroutine  (VB01AD)  from  the  Harwell  Subroutine 
Library(ref .8) .  In  order  to  remove  bias  introduced  by  the  concentration  of 
data  points  in  the  upper  layers  of  the  ocean,  it  was  found  necessary  to  weight 
the  data  samples  according  to  depth  so  that  the  correct  behaviour  was 
approached  in  the  deep  ocean.  Each  data  point  was  multiplied  by  a  weighting 
factor  determined  as 


W  =  1000(z/5000)n 


(4) 


where  the  exponent  n  was  chosen  by  eye  to  give  the  desired  behaviour  for  each 
area.  Values  of  the  parameter  so  computed,  together  with  the  weighting 
exponents,  are  listed  in  Appendix  III.  The  resulting  density-depth 
relationships  are  shown  by  the  solid  curves  in  figure  4. 

3.  THE  PROPERTIES  OF  A  COLUMN  OF  SEAWATER 

Given  that  the  temperature  and  salinity  profiles  for  a  column  of  seawater  are 
known,  then  by  invoking  suitable  relationships  it  is  possible  to  estimate  a 
large  number  of  the  physical  properties  of  the  seawater  column.  A  number  of 
such  relationships  are  introduced  and  discussed  in  the  present  section.  The 
intention  has  been  to  select  relationships  which  offer  a  high  degree  of 
precision  and  accuracy  and,  wherever  possible,  have  wide  acceptance  among  the 
scientific  community.  The  discussion  presupposes  the  incorporation  of  the 
techniques  and  relationships  described  into  a  digital  computational  scheme. 
In  many  places  therefore,  a  degree  of  sophistication  and  complexity  has  been 
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adopted  which  would  he  unacceptable  to  a  user  without  computing  facilities. 
Tables  of  the  properties  of  seawater  at  selected  values  of  temperature, 
pressure  and  salinity,  computed  using  the  relationships  derived  below,  are 
given  in  Appendix  IV. 

3.1  Gravitational  acceleration 

The  gravitational  acceleration  at  sea-level  may  be  regarded  as  a  function 
of  geographical  latitude  only,  and  is  given  by  the  expression(ref . 13) : 


go  =  9.78049(1  +  0. 0052884s in2(J)  -  0 . 0000059sin22<|>)  (5) 

where  0  is  the  latitude  and  g0  is  given  in  units  of  m.s-2.  If  z  represents 
depth  (in  metres)  measured  positive  downwards  from  the  ocean  surface,  then 
a  correction  may  be  made  to  account  for  departure  from  the  ocean  surface: 


g  =  g0(l  +  2.26Xl0"1z) 


(6) 


3.2  The  equation  of  state  for  seawater 

The  equation  of  state  for  seawater,  which  expresses  specific  volume  as  a 
function  of  pressure,  salinity  and  temperature,  has  been  the  subject  of 
extensive  investigation  in  recent  years.  A  critical  review  of  past 
studies,  and  various  forms  of  the  equation  of  state,  is  contained  in 
reference  14.  At  the  present  time,  the  introduction  of  a  new  standard 
equation  of  state  is  imminent(ref . 15) ,  the  format  of  the  equation  of  state 
following  that  of  previous  expressions  developed  in  references  9  and  14. 
The  computation  of  specific  volume  given  temperature,  pressure  and  salinity 
may  be  seen  as  a  three-stage  process.  Firstly,  the  density  of  pure  water 
is  computed  from  the  specified  temperature.  The  perturbation  to  the 
density  introduced  by  the  presence  of  salinity  is  then  computed.  Chen  and 
Millero(ref . 14)  have  pointed  out  that  it  is  advantageous  from  the  point  of 
view  of  estuarine  and  groundwater  studies  to  consider  pure  water  as 
infinitely  diluted  seawater  in  this  manner,  rather  than  deriving  a  one- 
atmosphere  equation  of  state  for  seawater  having  validity  over  the  limited 
range  of  salinities  which  may  be  encountered  in  the  deep  ocean.  Finally, 
the  specific  volume  at  pressure  may  be  computed  from  the  density  at  one 
atmosphere.  In  the  new  equation  of  state,  the  second-order  bulk  secant 
modulus  equation  is  used  for  this  purpose.  That  is, 

K  =  PV°/(V°  -  VP)  =  K°  +  AP  +  BP2  (7) 


where  K  is  the  secant  bulk  modulus  at  applied  pressure  P  (bars),  K°  is  the 
secant  bulk  modulus  at  one  atmosphere  (P=0),  V°  is  the  specific  volume 
(cm3.g_1)  at  atmospheric  pressure,  and  VP  is  the  specific  volume  at 
pressure  P.  A  and  B  are  adjustable  parameters.  As  of  April  1980  (the 
date  of  publication  of  reference  15),  the  one-atmosphere  equation  of  state 
remained  to  be  finalised.* 


*  Note  added  in  press.  Since  completion  of  the  work  reported  herein,  a 
new  international  standard  one-atmo sphere  equation  of  state  for  seawater 
has  appeared  in  the  literature.  Details  will  be  found  in  reference  26. 
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Recommendations  for  the  computation  of  K°  ,  A  and  B  in  equation  (7)  had  j 

previously  been  proposed  by  the  .Joint  Panel  on  Oceanographic  Tables  and 

Standards,  these  being  published  in  reference  lb.  This  source  would 

appear  however,  to  contain  typographical  errors,  since  computations  with 

the  expressions  contained  therein  failed  to  reproduce  the  results  of 

references  9  and  14  to  within  the  accuracy  stated  in  reference  16.  As  an 

interim  measure  therefore,  pending  the  publication  of  the  new  standard 

equation  of  state  in  more  positive  form,  the  sound-speed  derived  equation 

of  state  of  Chen  and  Millero(ref . 14)  has  been  adapted.  There, 


K°  =  K°  +  aS  +  bS1  *s  (7(a)) 

w 


A  =  Aw  +  cS  +  dS1  ’5  (7(b)) 


B  =  B  +  eS  (7(c)) 

w  v  ' 

where  the  salinity  S  is  in  parts/thousand,  and  the  subscript  w  refers  to 
pure  water.  The  parameters  a,  b,  c,  d  and  e  are  temperature  dependent  and 
given  by, 

a  =  53.751  -  0.4607T  +  7.  OSO.xlO-3'!'2  -  5.107xl0-sT3  (8) 


b  =  0.2322  -  4 . 838xl0-3  T  (9) 


c  =  4 . 69 2x1 0-3  -  8. 387xl0-5  T  +  4. 68x10-’ T2 


(10) 


d  =  -1.332x10 


(11) 


e  =  -1.4 12xl0-6  +  9.006xl0-8  T  -  1.551xlO-9TJ 


(12) 


where  T  is  the  temperature  in  degrees  centigrade.  The  equation  of  state 
for  pure  water  derives  from  the  work  of  Chen,  Fine  and  Millero(ref. 17) . 
The  pure  water  density  at  specified  temperature  and  atmospheric  pressure, 
and  associated  pressure  correction  parameters  are  computed  using. 
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p°  =  0.9998395  +  6.79l4xiO-5T  -  9 . 0894xl0-6T2 
w 

+  1 ,0171X10-7T3  -  1.2846x]0-9T4 

+  1 . 1592X10-1  ‘T5  -  5 . 0125x10-  14T6  (13) 

K°  =  19652.17  +  148.183T  -  2.29995T2  -  0.012810T3 
w 

-  4.91564xio-sT4  +  1 . 03553xi0-7T5  (14) 

A  =  3.26138  +  5.223xl0-4T  +  1.324xl0-4T2 
w 

-  7.655X10-7T3  +  8.584xiO-10T4  (15) 


B  =  7.2061X10-5  -  5.8948xi0-6T  +  8.699xi0-8T2 
w 

-  1.010X10-9T3  +  4.322X10-12T4  (16) 


where  p°  is  the  one  atmosphere  density  for  pure  water  (cm3g  *),  defining 
the  specific  volume, 


V°  =  1/Pj  07) 

The  perturbation  to  be  applied  to  p°  to  evaluate  the  density  of  seawater  at 
one  atmosphere  pressure  (p°)  derives  from  the  work  of  Mi  Hero,  Gonzalez  and 
Ward(ref . 18) .  w 
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p°  -  p9  =  (8 . 259 1 7*10-4  -  4 .  4490*  10-6T  +  1 . 0485X10-7T2 
-  1 .2580xi0-9T3  +  3.315X10- 12T4 )S  + 
(-6.33761X10-6  +  2.844lxlO-7T  -  1 .687lxl0-sT2 


+  2 . 83258x10- 10T3)S 1  ,s  + 

(5.4705X10-7  -  1.97975xiO-®T  +  1 .6641xlO-9T2 
-  3. 1203X10-1 *T3)S2  (18) 


and 


V°  =  l/p° 


(19) 


The  equation  of  state  of  seawater  contained  in  equations  (7  to  19)  may  be 
used  to  compute  specific  volume  to  an  accuracy  of  (±10xl0-6  cn^.g-1)  over 
the  range  of  0  to  40  parts/thousand  in  salinity,  0  to  40°C  and  0  to  1000 
bars.  In  addition,  the  equation  of  state  can  be  used  to  determine 
isothermal  compressibility  to  an  accuracy  of  ±0.03xl0-6  bar-1  and  thermal 
expansibility  to  an  accuracy  of  ±2xl0-6  degrees-1  over  the  same  range  of 
salinity,  temperature  and  pressure. 

The  isothermal  compressibility,  which  finds  application  in  the 
determination  of  seawater  density  from  in  situ  measurements,  and  in  the 
determination  of  sound  speed(ref . 19) ,  is  defined  as. 


p  e  — L(avP/3P) 
V 


(20) 


and  may  be  determined  from  equation  (7) : 


V°(K°  -  BP2) 
VP(K°  +  AP  +  BP2)2 


(21) 


Similarly,  the  expansibility  of  seawater 
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a  a  -i-OVP/3T) 

V 

can  be  determined  as 

POV/3T)0 

a  =  -^(3V/3T)° - p - 

V  V  (K°  +  AP  +  BP2) 

(3K°/3T)d  +  P(3A/3T)_  +  P2(3B/3T)D 

+  PV° - - - - - — ^ - — 

V  (K°  +  AP  +  BP2)2 

wherein  the  derivatives, 

(3V/3T)0  =  -(V°)2(3p/3T)° 

(3p/3t)°  =  6.7914X10-5  -  1 .81788*10-5T 

+  3.0513xl0-7T2  -  5.1384xl0-9T3 
+  5 . 796x10- 1 1T4  -  3 . 0075X10- 13T5 
+  (-4.4490x1 0-6  +  2.097x1 Q-7T 

-  3.774xi0-9T2  +  1.326xlO-11T3)S 
+  (2.8441X10-7  -  3.3742xlO-8T 

+  8 . 49774x10- j0T2)S1 ' 5  +  (-1.97975X10-8 


(22) 


(23) 


(23(a)) 


+  3.3282xi0-9T  -  9.3609xi0-nT2)S2 


(23(b)) 
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(3K°/3T)  =  148.183  -  4.5999T  +  3.843*10-2T2 

-  1.966256xi0-4T3  +  5 . 17?65xi0-7T4 
+  (-0.4607  +  1 . 406*10-2T  -  1  ..r>32lxiO-4T2)S 
+  (-4.838xlO-3)S1 ' 5 

(23(c)) 

(3A/3T)p  =  5 . 223* 10-4  +  2.648xi0-4T  -  2 . 2965*10-6T2 

+  3.4336xi0-9T3  +  (-8.387X10-5  +  9.36xiO-7T)S  (23(d)) 


(3B/3T)p  =  -5.8948X10-6  +  1.7398xl0-7T  -  3.03xiO-9T2 

+  1 .7288x10-  1  *T3  +  (9.006x10-*  -  3 . 102xi0-9T)S  (23(e)) 


3.3  Pressure 

Employing  the  customary  units  as  detailed  in  Sections  3.1  and  3.2  above, 
the  pressure  at  depth  Z  (metres),  measured  positive  downwards,  may  be  found 
by  the  relationship, 


P  = 


Z 

1 0-2J  gpdz 
0 


(24) 


Using  the  one-atmosphere  density  profile  computed  using  equations  (13) 
and  (18)  as  a  first  approximation  to  p,  a  cubic  spline  is  fitted  through 
the  recorded  values  of  the  product  g(z).p(z)  using  the  Harwell  Subroutine 
TB04AD(ref .8) .  The  cubic  spline  is  integrated  between  the  surface  and  the 
appropriate  depths  using  Harwell  Subroutine  QG02AD  to  give  a  first 
approximation  to  the  pressure  profile  using  equation  (24).  Most  recent 
estimates  of  the  density  and  pressure  profiles  are  then  used  to  iteratively 
update  one  another. 

3.4  The  speed  of  sound 

The  speed  of  sound  U,  in  a  compressible  medium  is  defined  by  the  equation, 


1/U2  =  (3p/3P)0  =  - 1 /V2 (3V/3P)0 


(25) 


-  1 1 
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where  0  is  the  entropy.  By  expanding  the  right-hand  side  of 

equation  (25),  this  may  be  rewritten  as(ref.20), 

1/U2  =  -1/V2[(9V/3P)t  +  T' OV/3T)2/Cp]  (2b) 

where,  in  this  context,  T'  is  understood  to  refer  to  absolute  temperature. 
0^  is  the  specific  heat  at  constant  pressure,  and  p  =  -1/V  (9V/9P)^,  is  the 

iso^ermal  compressibility  previously  introduced  (Section  5.2). 
Equation  (26)  then,  comprises  a  unique  relationship  between  sound  speed, 
the  equation  of  state  and  the  specific  heat,  and  forms  the  basis  for 
methods  of  deriving  the  equation  of  state  of  water  based  on  sound  speed 
measurements.  In  equation  (26),  the  first  term  is  the  dominant  term,  the 
thermodynamic  term  forming  a  second-order  correction(ref .20) .  In  this  and 
the  next  section,  expressions  for  the  speed  of  sound  in  seawater,  and  the 
specific  heat  are  presented. 

The  form  of  the  relationship  equation  (26)  establishes  the  desirability  of 
achieving  a  level  of  consistency  between  the  equation  of  state  of  seawater 
used  in  a  particular  study,  and  the  sound  speed  equation.  For  this  reason, 
the  sound  speed  equation  of  Chen  and  Millero(ref .20)  has  been  adopted  in 
the  present  study,  this  forming  the  basis  upon  which  their  equation  of 
state(ref . 14)  was  derived.  The  sound  speed  equation  consists  of 

relationships  for  the  speed  of  sound  in  pure  water,  and  for  the 
perturbations  introduced  by  the  presence  of  salinity.  The  composite  system 
of  equations  is  valid,  as  is  the  equation  of  state,  over  a  range  of  0  to  40 
parts/thousand  in  salinity,  0  to  40°C,  and  0  to  1  000  bars  in  pressure. 
In  the  following,  sound  velocities  are  measured  in  m.s-1,  other  variables 
being  described  in  terms  of  their  conventional  units. 

For  pure  water,  the  speed  of  sound  is  expressed  as  a  function  of 
temperature  and  pressure, 


U  =  1402.288  +  5.0371  IT  -  5 . 80852X 10-2T2 

w 

+  3 . 3420* 10-4T3  -  1 .47800xi0-6T4  +  3 . I464xi0-9T5  + 
(0.153563  +  6.8982X10-4T  -  8. 1788xiO-6T2 
+  1 .3621X10-7T3  -  6. 1185xl0-10T4)P  + 

(3.1260X10-5  -  1.7107xl0-6T  +  2 .5974xi0-8T2 
-  2. 5335x10- 10T3  +  ] ,0405X10-12T4)P2  + 


(-9.7729X10-9  +  3 . 8504x10- 1 °T  -  2.3643xi0-12T2)P3 


(27) 
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At  atmospheric  pressure,  the  differential  between  the  speed  of  sound  in 
pure  water,  and  that  in  seawater  is  given  by, 


II  -  U°  =  (1.189  -  1 . 262*  10-2T  +  7 . I64x  10-5T2 
w  w 


+  2.006xl0-6T3  -  3.21X10-8T4)S  + 


(-1.922X10-2  -  4.42X10-5T)S1 ' 5  +  1.727xiO-3S2  (28) 


From  the  above,  the  speed  of  sound  in  seawater  under  pressure  may  be  found 
using  the  relationship, 


(UP  -  UP)  -  (U°  -  U°)  =  A  S  +  B  S1’5  +  C  S2 
w  w  s  s  s 


(29) 


where 


A  =  (9.4742X10-5  -  1.2580xl0-sT  -  6 . 4885xio-8T2 
s 

+  1 . 0507X 10-8T3  -  2 .0122X10- 10T4)P  + 

(-3.9064X10-7  +  9.1041X10-9T 

-  1 ,6002X10-10T2  +  7.988xi0-12T3)P2  + 

( 1 .  100x10- 10  +  6 . 649x10- 12T 

-  3. 389x10- 13T2)P3  (29(a)) 


Bg  =  (7.3637X10-5  +  1 . 7945xi0-7T)P 


(29(b)) 


C  =  -7.9836xiO-6P 
s 


(29(c)) 


Finally,  it  should  be  mentioned  that  the  work  of  Kroebel(ref . 14,  Annex  4) 
indicates  that  the  presence  of  suspended  particles  may  be  a  signficant 
additional  factor  in  the  determination  of  sound  speed  under  operational 
conditions . 
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3.5  Specific  heat 

The  specific  heat  of  seawater  is  of  importance,  not  only  on  account  of  its 
connection  with  the  equation  of  state,  and  the  sound  speed  equation,  but  in 
its  own  right,  as  a  vital  component  in  the  thermodynamic  make-up  of  the 
oceanic  environment.  As  such,  it  is  especially  important  in  the 
computation  of  adiabatic  lapse  rates(refs . 16 ,20) .  While  a  number  of 
investigations  of  the  specific  heat  of  seawater  at  atmospheric  pressure 
have  been  performed(refs.22,23) ,  there  are  in  existence  no  measurements  of 
the  specific  heat  of  seawater  at  elevated  pressures(ref . 16) .  If  however, 
the  heat  capacity  can  be  evaluated  at  atmospheric  pressure  for  specified 
temperature  and  salinity,  then  the  specific  heat  at  elevated  pressure  can 
be  computed  using  the  theoretical  relationship(ref .24) , 


P 

>  =  Cp  ‘  T’/  (32V/3T2)pdP  (30) 


where  in  this  equation,  T'  is  understood  to  be  the  absolute  temperature. 

Millero,  Perron  and  Desnoyers(ref .23)  give  a  relationship  for  the  specific 
heat  of  seawater  at  atmospheric  pressure, 


C°  =  C°  +  A  S  +  B  S1-5  (31) 

p  pw  cp  cp 

where  C°  is  the  specific  heat  (j.g-1  .deg-1  )  of  seawater,  and  C°  is  the 
p  P" 

specific  heat  of  pure  water  at  the  same  temperature,  which  may  be  found 
from, 

C°  =  4.2174  -  3 . 720283*10-3T  +  1 . 412855*10-4T2 
pw 

-  2.654387*10-6T3  +  2.093236xl0-8T4  (31(a)) 


and  the  temperature  dependent 


parameters  A  and  B 
r  cp  cp 


are  given  by, 


A  =  -7.644X10-3  +  1.0727xi0-4T  -  1.38xi0-6T2 
cp 


(31(b)) 


B  =  1.77X10-4  -  4.08X10-6T  +  5.35xl0-8T2  (31(c)) 

cp 

The  pressure-dependent  term  of  equation  (30)  may  be  derived  from  the 
equation  of  state  of  seawater  (Section  3.2).  That  is,  if  we  let, 


P 

1 2  =  S  02V/3T2) pdP, 
0 
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then,  equation  (30)  may  be  rewritten  as, 

Cp  =  C£  -  0. 1(T  +  273. 15 ) I 2  (32) 

where , 

h  =  PCa^/at2)0  +  A,  I,,  +  B,  I21  +  feljj  +  B3Ij3 

+  C1I3X  +  C2I32  +  C3I33  +  C4I34  +  C5I35  (33) 

If  we  let 

X  =  K°  +  AP  +  BP2  (34(a)) 

and 

q  =  4K°B  -  A2  (34(b)) 

then  the  integrals  represented  by  the  I  terms  above  may  be  evaluated  as 
follows : 

P 

1 10  =  J  =  2[tan-1{(2BP+A)/Vq}  -  tan-1 {A/Vq} ]/Vq»  (q>0) 

0 

=  log[  (2BP  +  A  -  ^q)(A  +  Vrq)/{(2BP  +  A  +  ^q)(A  -  ^q)}]/^q, 

(q<0)  (34(c)) 

P 

In  =  J  -|-dP  =  -2^” ln(X/K°)  -  -A^  (34(d)) 


I20  =  /  41-  =  (2BP+A)/ qX  -  A/qK°  +  2BI,0/q  (34(e)) 

0  * 
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hs  =/  Ja  dp  =  yl23  -  g-  I33  -JI34  (34 (n) ) 

The  coefficients  in  equation  (34)  are 

A,  =  -^V/aT3)0  (35(a)) 

Bi  =  2(3V/3T)°  (3K°/3T)p  +  V°  (92  K°/aT2)p  (35(b)) 

Bj  =  2(3V/3T)°  (3A/3T)p  +  V°  (32A/3f!  )p  (35(c)) 

B3  =  2(3V/3T)°  (3B/3T)p  *  V°  (32  B/3T2 )p  (35(d)) 

C,  =  -2V°  (3K°  /3'  "'2  (35(e)) 

T 

C2  =  -4V°(3K°/3T)p(3A/3T)p  (35(f)) 

C3  =  -4V°  (3K° /3T) p(3B/3T) p  -  2V°  (3A/3T)3p  (35(g)) 

C«  =  -4V°  (3A/3T)p(3B/3T)p  (35(h)) 

C*  =  -2V°  (3b/3T)p  (35(i)) 


In  the  above,  V°  and  the  parameter  ,  A  and  B  are  given  by  equations 
(7)  to  (19).  The  first  derivatives  are  computed  according  to 
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equations  (23(a)) to  23(d)),  and  the  second  derivatives  may  be  found  using, 

(32V/3T2)°  =  -(V°)2(32p/3T2)° 

+  2(V°)3[( 3p/3T)°J2  (36(a)) 

(32p/3T2)0  =  -1.81788X10-5  +  6.1026xio-7T  -  1 .54152X10-8T2 
+  2,3184xlO-10T3  -  1 .50375xi0-12T4 
+  (2.097X10-7  -  7.548xiO-9T  +  3.978xiO-11T2)S 
+  (-3.3742X10-8  +  1.699548xl0-9T)S1'5 
+  (3.3282X10-9  -  1.87218xlO-10T)S2  (36(b)) 

(32K°/3T2)p  =  -4.5999  +  7.686xlO-2T  -  5 .898768xi0-4T2 

+  2.07106xiO-6T3  +  (1.406X10-2  -  3 . 0642xi0-4T)S  (36(c)) 

(32A/3T2)p  =  2.648X10-4  -  4.593xiO-6T 

+  1.03008xi0-*T2  +  9.36xiO-7S  (36(d)) 

(32B/3T2)p  =  1.7398X10-7  -  6.06xi0-9T 

+  5 . 1864x10- 1 ‘T2  -  3.102X10-9S  (36(e)) 

Computer  evaluation  of  I2  using  equations  (34)  and  (36)  becomes  inaccurate 
if  |B|  becomes  too  small.  As  a  check  on  the  accuracy  of  the  evaluation  of 
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la*  this  integral  was  evaluated  numerically.  That  is,  the  expression 

(32V/3T2)p  =  (32V/3T2)0  -  P(3V/3T)°/X  +  PV0Xi/X2 

+  PV°X2/X2  -  2PV°Xj/X3  (37(a)) 


where 


X2  =  (3K°/3T)p  +  P(3A/3T)p  +  P2(3B/3T)p  (37(b)) 


and 


X3  =  (32K°/3T2)p  +  P(32A/3T2)p  +  P2(32B/3T2)p  (37(c)) 


was  integrated  between  the  limits  of  0  and  P  using  an  adaptive  quadrature 
routine  (Harwell  subroutine  QA05AD) .  Using  double  precision,  the  results 
computed  using  the  two  approaches  were  found  to  agree  to  five  figures 
provided  |B|^5*10-6.  The  approach  adopted  for  the  computation  of  specific 
heat  is  a  hybrid  one;  when  |B|g5><10-6,  I2  is  evaluated  analytically,  while 
for  | B | <5X 1 0— 6  the  numerical  scheme  described  above  is  used. 

Even  assuming  that  I2  can  be  determined  precisely,  in  the  absence  of 
measurements  of  specific  heat  of  seawater  at  elevated  pressures,  and  as  a 
result  of  the  shortage  of  measurements  at  one  atmosphere  also,  the 
precision  with  which  equations  (29)  and  (30)  may  be  used  to  estimate  the 
specific  heat  of  seawater  remains  in  some  doubt.  Millero(ref . 16,  p  14) 
however,  comments  on  accuracy  of  the  one-atmosphere  relationship  for 
specific  heat  (equation  (30)),  and  feels  that  it  is  sufficiently  accurate 
for  the  computation  of  lapse  rates,  but  possibly  not  for  the  definition  of 
entropy.  He  states  furthermore,  that  measurements  of  specific  heat  of 
pure  water  at  elevated  pressures  do  confirm  the  validity  of  using  the 
thermodynamic  relationship  (equation  (29))  to  predict  values  at  elevated 
pressures . 

3.6  Dynamic  height  anomaly 

Dynamic  height  anomaly  provides  a  measure  of  the  mass  distribution  in  an 
ocean  water  column,  its  theoretical  significance  and  use  in  the  computation 
of  geostrophic  currents  being  described  in  most  texts  on  oceanography 
(reference  13,  chap  8;  reference  14,  chap  7).  Dynamic  height  anomaly  is 
defined  as 


P 

AD  =  f  6dP 
Po 


(38) 


11) 


where  l’0  is  the  pressure  ;it  the  ocean  surface  (Po=0),  and  I’  is  the  pressure 
at  depth.  6  is  the  specific  volume  anomaly, 

6  =  V(S.T.P)  -  V (35,0, P)  (39) 


which  provides  a  measure  of  the  departure  of  the  specific  volume  of  the 
ocean  from  that  of  a  "standard  ocean"  having  salinity,  S=3S  parts/thousand, 
temperature,  t=0  (',  and  pressure  P.  In  the  appended  computational  package 
the  integration  of  equation  (39)  is  performed  by  fitting  a  cubic  spline  to 
5=f(P)  using  the  Harwell  library  subroutine  TB04BD,  and  integrating  the 
resulting  piecewise-cubic  polynomial  using  the  Harwell  subroutine  QG02AD. 

One  of  the  fundamental  problems  associated  with  the  computation  of  ocean 
currents  using  dynamic  height  anomaly  is  that  of  establishing  the  absolute 
current  values;  the  dynamic  method  provides  relative  currents  at  different 
depths  only.  As  an  extension  of  the  dynamic  method*  one  may  seek  a 
regression  of  on  T _  j / 2 2  w^cre  z *  an  assumed  level  of  no  motion, 

and  zl  is  some  depth  less  than  z2;  *s  t*ie  dynamic  height  anomaly 

between  zl  and  z2,  being  the  temperature  difference  over  this  depth 

increment.  Now,  if  z2  is  a  level  of  no  motion,  it  follows  that  within  a 
certain  geographical  area  t  ^  will  be  approximately  constant,  and  we  may 

alternatively  seek  a  regression  of  D  .  on  T  . .  Provided  z2-zl  is 

Zl/ZZ  Zl 

sufficiently  small,  it  is  reasonable  to  expect  that  D  _  will  correlate 

Z 1  /  Z  2. 

closely  with  T  Historically,  for  the  East  Australian  Current  area,  z2 

is  set  at  1300  m  while,  in  the  interests  of  accuracy,  zl  should  not  be 
taken  as  less  than  450  m.  An  attempt  to  measure  the  regression  of 
D450/ 1300  on  T450  was  made  during  the  present  study,  but  the  available  data 

base  proved  too  small  to  allow  reliable  relationships  to  be  computed. 
Pearce  and  Hamon(ref.7)  have,  however,  provided  regression  curves  giving 
partial  coverage  of  the  area  under  consideration.  An  example  of  the  use 
of  such  relationships  is  given  by  reference  25. 


4.  COMPUTER  PROGRAMME  ORGANISATION 

A  computer  programme  has  been  written  to  compute  and  plot  the  parameters 
discussed  in  the  preceding  section  for  a  series  of  XBT  profiles.  The 
programme  is  of  a  modular  form  comprising  a  number  of  separate  routines  which 
be  used  independently  if  so  desired.  The  progranme  has  been  written  to  run 
on  the  IBM  370/3033  computer  at  DRCS  and  provide  graphical  output  on  a  CALC0MP 
12-inch  plotter.  ANSI  standard  l-'ORTRAN  IV  has  been  used  throughout,  but  the 
package  calls  on  subroutines  from  two  external  libraries.  Subroutines  from 
the  Harwell  Subroutine  Library  have  been  used  to  fit  cubic  splines  to  data 
points  and  to  perform  interpolation  and  quadrature  using  the  fitted  curves. 


*  As  pointed  out  by  a 
S.  Godfrey  (CSIR0) . 


referee. 


this  method  was  originally  suggested  by 


A. 
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and  for  numerical  quadrature.  If  this  library  is  unavailable  it  is  a 
straightforward  matter  to  replace  these  subroutines  with  suitable  equivalents. 
To  facilitate  this  procedure,  details  of  the  calls  to  the  subroutines  used  are 
contained  in  Appendix  IX.  In  addition,  a  number  of  CALCOMP  plotter 
subroutines  are  used,  and  since  these  may  not  be  available  for  use  with 
another  installation,  details  of  the  appropriate  calls  are  given  in 
Appendix  X.  An  overview  of  the  programme  appears  in  figure  5,  and  a  programme 
listing  presented  in  Appendix  V. 

For  each  XBT  profile  processed  the  programme  produces  a  title  block  containing 
environmental  data,  followed  by  a  listing  in  tabular  form  of  the  computed 
oceanographic  parameters  associated  with  the  XBT  profile.  A  sample  output  is 
presented  in  Appendix  VI.  In  addition  the  programme  optionally  produces 
graphical  output  of  two  types.  The  first  type  of  plot  (figure  6;  hereafter 
referred  to  as  a  'fourplot')  shows  temperature/depth,  salinity/depth, 
sound  speed/depth  and  sigma-t/depth  profiles  for  a  specified  XBT  drop.  In 
the  second  type  of  plot  (figure  7;  hereafter  referred  to  as  a  'longplot'), 
temperature/depth  or  sound  speed/depth  profiles  for  a  sequence  of  XBT's  are 
plotted. 

4.1  Plotting  options 

The  selection  of  the  plotting  options  is  made  by  means  of  a  list-directed 
record  input  on  logical  unit  5.  The  various  options  are  specified  by  four 
integer  variables  as  follows: 


First  integer 
second  integer 

third  integer 
fourth  integer 


1  to  produce  fourplots, 

1  for  shallow  XBT  (T-4)  type, 

2  for  deep  XBT  (T-7)  type, 

1  if  any  continuous  longplots  are  required, 

1  if  a  temperature  longplot  is  required, 

2  if  a  sound-speed  longplot  is  required, 

3  if  both  temperature  and  sound-speed  logplots 
are  required. 


Hence,  a  plot  option  record  of  1,1, 1,1  would  cause  a  fourplot  and  a 
longplot  of  temperature  versus  depth  to  be  plotted  with  the  depth  axis 
scaled  for  T-4  type  XBT's.  A  blank  record  read  on  logical  unit  5  will 
suppress  plotter  output.  If  any  longplots  are  required,  a  temporary  work 
file  must  be  allocated  to  logical  unit  3. 


4.2  Temperature/Depth  data  input 


The  programme  is  written  to  accept  data  from  a  card  image  file  input  on 
logical  unit  2.  The  following  data  format  has  been  in  use  within  Marine 
Sciences  Composite  (MSC)  at  DRCS  since  1975.  All  XBT  data  gathered  by  M SC 
and  stored  on  the  IBM  370/3033  have  been  written  in  this  format,  and  have 
been  validated  and  calibrated.  A  eard  image  file  to  process  a  series  of 
XBT  profiles  using  the  appended  programme  should  be  written  in  the 
following  format: 


J 


I 


WSKI  ’‘i-TM 


KKCOKI)  I:  Contains  the  number  ol  XBT  profiles  to  tollow  in  110  iorm.il, 
left  j us t i f ied . 

RECORD  2:  A  36  character  record  (9A4),  left  justified,  containing  the 
cruise  title. 

RECORD  3: 


Element 

Starting 

column 

Format 

Day 

1 

A2 

Month 

4 

A2 

Year 

7 

A2 

T  ime 

12 

A2 

Longitude,  °E 

20 

F7.2 

Latitude,  °S 

38 

F7.2 

RECORD  4: 


Element 

Starting 

column 

Format 

Sonic  Depth,  M. 

2 

A5 

Wet  Temperature,  °C 

17 

A4 

Dry  Temperature,  °C 

26 

A4 

Wind  Direction 

36 

A3 

Wind  Speed,  knots 

39 

A3 

Barometric  Pressure,  mb. 

49 

A6 

RECORD  5: 


Element 

Starting 

Format 

column 

Profile  Number 

1 

13 

No.  of  Depth/Temperature 

readings  (L) 

4 

15 

RECORDS  6  to  end  of  profile: 

Contain  a  sequence  of  L  pairs  of  depth  and  temperature  readings  in  the 
order  D(l) ,T(1) ,D(2) ,T(2) , . . . ,D(l ) ,T(L)  packed  in  12F6.2  format  in  80 
character  records. 

Records  2  to  6+  are  repeated  for  each  profile  in  the  batch. 

A  sample  input  file  of  this  type  is  included  in  Appendix  VI 1. 

The  above  input  format  is  applicable  only  to  data  originating  from  Marine 
Sciences  Composite;  users  wishing  to  process  data  from  other  sources  must 
make  the  appropriate  changes  to  the  programme  input.  Details  of  the 
standard  format  adopted  by  the  Australian  Oceanographic  Data  Centre  for 
storage  of  XBT  data  within  their  current  data  base  are  given  in 
Appendix  VIII.  Note  that  the  AODC  format  does  not  store  certain 

components  of  the  environmental  data  which  are  stored  by  MSC;  in  order  to 
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itli  list-  the  output  format  used  liy  the-  existing  programme  it  is  only 
necessary  to  extract  cruise  title  and  positional  information  together  with 
depth  and  temperature  figures  from  the  input  records.  Note  also  that 
temperature  and  depth  data  supplied  by  AO DC  will  also  require  correction 
for  offset;  the  necessary  calibration  information  is  supplied  by  the 
parameters  CALDEP  and  CALTEM  in  the  input  records. 


i 
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APPENDIX  I  NANSKN  CASTS  USED  Tt»  DERIVE  T-S  CURVES 


CAST 

SOURCE  CRUISE 

CAST 

SEASON  AREA 

DATE 

LATITUDE 

LONGITUDE 

DEPTH 

NO. 

I DENT. 

NO. 

NO. 

CODE 

DEGREES 

DEGREES 

M. 

PTS 

1 

RANRL 

S4/3 

36 

WINTER 

A1 

19/07/68 

18.00 

S 

152.35 

E 

1756 

23 

■> 

CSIRO 

G 1/62 

11 

SUMMER 

A1 

18/01/62 

18.62 

S 

155.97 

E 

1981 

18 

3 

CSIRO 

Gl/62 

12 

SUMMER 

A1 

19/01/62 

18.33 

S 

157.72 

E 

3450 

20 

4 

CSIRO 

Gl/62 

14 

SUMMER 

A1 

20/01/62 

22.97 

S 

156.95 

E 

1983 

15 

5 

CSIRO 

Gl/62 

15 

SUMMER 

A1 

20/01/62 

22.90 

S 

158.70 

E 

2465 

16 

6 

CSIRO 

Gl/62 

24 

SUMMER 

A1 

27/01/62 

24.67 

S 

158.33 

E 

2428 

16 

7 

CSIRO 

Gl/62 

26 

SUMMER 

A1 

28/02/62 

24.70 

S 

155.37 

E 

1956 

14 

8 

RANRL 

S4/3 

13 

WINTER 

A1 

13/07/68 

25.00 

S 

154.00 

E 

1928 

24 

9 

RANRL 

S4/3 

17 

WINTER 

A1 

14/07/68 

25.12 

S 

157.00 

E 

1871 

23 

10 

RANRL 

S4/3 

20 

WINTER 

A1 

15/07/68 

22.92 

S 

154.95 

E 

1977 

21 

1 1 

RANRL 

S4/3 

25 

WINTER 

A1 

16/07/68 

20.75 

s 

153.50 

E 

1564 

22 

12 

RANRL 

S4/3 

30 

WINTER 

A 1 

17/07/68 

20.67 

s 

157.03 

E 

1796 

23 

13 

CSIRO 

G3/63 

124 

WINTER 

A2 

16/08/63 

24.07 

S 

169.30 

E 

1500 

13 

14 

CSIRO 

Gl/60 

74 

SUMMER 

A2 

22/02/60 

18.17 

S 

167.43 

E 

4410 

20 

15 

CSIRO 

Gl/60 

66 

SUMMER 

A2 

22/02/60 

15.00 

s 

163.01 

E 

4280 

20 

16 

CSIRO 

Gl/60 

69 

SUMMER 

A2 

22/02/60 

17.20 

s 

163.62 

E 

3307 

18 

17 

CSIRO 

Gl/60 

78 

SUMMER 

A2 

24/02/60 

20.52 

s 

169.02 

E 

2421 

16 

18 

CSIRO 

Gl/60 

80 

SUMMER 

A2 

25/02/60 

21.98 

s 

167.63 

E 

1775 

15 

19 

CSIRO 

G2/60 

112 

AUTUMN 

A2 

24/03/60 

23.32 

s 

165.05 

E 

3312 

19 

20 

CSIRO 

G2/60 

114 

AUTUMN 

A2 

24/03/60 

21.53 

s 

163.72 

E 

3410 

19 

21 

CSIRO 

G2/60 

119 

AUTUMN 

A2 

25/03/60 

17.92 

s 

161.17 

E 

2875 

19 

22 

CSIRO 

G2/60 

121 

AUTUMN 

A2 

26/03/60 

15.93 

s 

160.80 

E 

4193 

22 

23 

CSIRO 

Gl/62 

20 

SUMMER 

A2 

26/01/62 

24.62 

s 

166.42 

E 

2743 

17 

24 

CSIRO 

Gl/62 

23 

SUMMER 

A2 

27/01/62 

24.55 

s 

160.37 

E 

1700 

14 

25 

CSIRO 

Gl/60 

45 

SUMMER 

A3 

17/02/60 

17.55 

s 

176.12 

E 

2500 

17 

26 

CSIRO 

G3/63 

120 

WINTER 

A3 

14/08/63 

19.27 

s 

176.60 

E 

1968 

15 

27 

CSIRO 

G3/63 

121 

WINTER 

A3 

15/08/63 

20.38 

s 

174.85 

E 

2320 

14 

28 

CSIRO 

G3/63 

122 

WINTER 

A3 

15/08/63 

21.70 

s 

172.85 

E 

1972 

14 

29 

CSIRO 

G3/63 

123 

WINTER 

A3 

16/08/63 

22.88 

s 

171.07 

E 

2000 

14 

30 

CSIRO 

G3/63 

124 

WINTER 

A3/A2 

16/08/63 

24.07 

s 

169.30 

E 

1500 

13 

31 

CSIRO 

Gl/60 

38 

SUMMER 

A3 

10/02/60 

24.07 

s 

175.52 

E 

4273 

20 

32 

CSIRO 

Gl/60 

40 

SUMMER 

A3 

11/02/60 

22.75 

s 

175.97 

E 

3855 

14 

33 

CSIRO 

Gl/60 

41 

SUMMER 

A3 

11/02/60 

21.32 

s 

176.05 

E 

2000 

15 

34 

CSIRO 

Gl/60 

42 

SUMMER 

A3 

11/02/60 

20.78 

s 

176.33 

E 

2785 

17 

35 

CSIRO 

Gl/60 

44 

SUMMER 

A3 

12/02/62 

19.13 

s 

176.80 

E 

2977 

17 

36 

CSIRO 

Gl/60 

48 

SUMMER 

A3 

17/02/60 

16.68 

s 

174.48 

E 

2500 

17 

37 

CSIRO 

Gl/60 

50 

SUMMER 

A3 

18/02/60 

16.78 

s 

172.90 

E 

2425 

17 

38 

CSIRO 

G6/64 

242 

SPRING 

B1 

18/09/64 

31.82 

s 

157.90 

E 

1500 

14 

39 

RANRL 

S4/3 

6 

WINTER 

B1 

11/07/68 

29.17 

s 

158.03 

E 

1914 

24 

40 

CSIRO 

Gl/62 

36 

SUMMER 

B1 

30/01/62 

27.33 

s 

156.58 

E 

3365 

17 

41 

CSIRO 

G6/64 

263 

SPRING 

B1 

22/09/64 

34.93 

s 

151.40 

E 

1500 

14 

42 

CSIRO 

Gl/60 

8 

SUMMER 

B1 

04/02/60 

33.88 

s 

159.88 

E 

1948 

17 

43 

CSIRO 

G3/64 

137 

AUTUMN 

B1 

21/03/64 

34.95 

s 

154.25 

E 

1424 

14 

44 

CSIRO 

G4/63 

141 

SPRING 

B1 

11/09/63 

33.03 

s 

155.15 

E 

1463 

14 

45 

CSIRO 

G2/60 

91 

SPRING 

B1 

19/03/60 

30.40 

s 

156.53 

E 

4419 

19 

46 

CSIRO 

Gl/64 

35 

SUMMER 

B1 

20/01/64 

31.03 

s 

153.50 

E 

1500 

14 

47 

CSIRO 

Gl/64 

32 

SUMMER 

B1 

20/01/64 

28.58 

s 

154.15 

E 

1500 

14 

48 

CSIRO 

Gl/62 

23 

SUMMER 

B1 

27/01/62 

24.55 

s 

160.37 

E 

7000 

14 

49 

RANRL 

S4/3 

12 

WINTER 

B1 

12/07/68 

27.17 

s 

154.00 

E 

1591 

22 

27 
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50 

CSIRO 

Gl/62 

42 

SUMMER 

B2 

02/02/62 

27.62 

S 

170.05 

E 

1880 

13 

51 

CSIRO 

G2/60 

100 

AUTUMN 

B2 

21/03/60 

30.37 

S 

165.47 

E 

3240 

38 

52 

CSIRO 

Gl/62 

23 

SUMMER 

B2 

27/01/62 

24.55 

S 

160.37 

E 

1700 

14 

53 

CSIRO 

Gl/62 

38 

SUMMER 

B2 

31/01/62 

27.50 

S 

161.20 

E 

1600 

14 

54 

CSIRO 

Gl/60 

8 

SUMMER 

B2 

04/02/60 

33.88 

S 

159.88 

E 

1948 

17 

55 

CSIRO 

Gl/60 

13 

SUMMER 

B2 

05/02/60 

33.73 

S 

165.65 

E 

2432 

17 

56 

CSIRO 

Gl/60 

23 

SUMMER 

B2 

07/02/60 

33.82 

s 

169.10 

E 

1635 

15 

57 

CSIRO 

G2/60 

95 

AUTUMN 

B2 

20/03/60 

30.37 

s 

160.78 

E 

1410 

16 

58 

CSIRO 

G2/60 

105 

AUTUMN 

B2 

22/03/60 

28.18 

S 

168.38 

E 

3404 

19 

59 

CSIRO 

G2/60 

108 

AUTUMN 

B2 

23/03/60 

26.07 

S 

166.63 

E 

3300 

19 

60 

CSIRO 

G3/63 

126 

WINTER 

B2 

17/08/63 

26.20 

S 

165.82 

E 

3125 

17 

61 

CSIRO 

G3/63 

127 

WINTER 

B2 

17/08/63 

27.30 

s 

164.08 

E 

2000 

14 

62 

CSIRO 

Gl/62 

42 

SUMMER 

B3 

02/02/62 

27.62 

s 

170.05 

E 

1880 

13 

63 

CSIRO 

Gl/60 

36 

SUMMER 

B3 

10/02/60 

26.35 

s 

175.17 

E 

3900 

19 

64 

CSIRO 

Gl/oO 

25 

SUMMER 

B3 

08/02/60 

33.63 

s 

171.43 

E 

1300 

13 

65 

CSIRO 

Gl/61 

31 

SUMMER 

B3 

03/02/61 

34.80 

s 

170.38 

E 

2049 

15 

66 

CSIRO 

Gl/60 

32 

SUMMER 

B3 

09/02/60 

29.20 

s 

174.23 

E 

2842 

17 

67 

CSIRO 

Gl/62 

46 

SUMMER 

B3 

03/02/62 

32.07 

s 

173.37 

E 

2000 

14 

68 

CSIRO 

Gl/62 

44 

SUMMER 

B3 

02/02/62 

29.43 

s 

172.33 

E 

2000 

14 

69 

CSIRO 

Gl/62 

43 

SUMMER 

B3 

02/02/62 

27.67 

s 

171.95 

E 

3245 

18 

70 

CSIRO 

G2/60 

103 

AUTUMN 

B3 

22/03/60 

30.38 

S 

168.73 

E 

3696 

19 

71 

CSIRO 

Gl/60 

34 

SUMMER 

B3 

09/02/60 

27.65 

s 

174.77 

E 

3400 

16 

72 

CSIRO 

Gl/60 

30 

SUMMER 

B3 

09/02/60 

30.73 

S 

173.88 

E 

2400 

15 

73 

CSIRO 

Gl/60 

28 

SUMMER 

B3 

08/02/60 

32.20 

s 

173.50 

E 

2475 

16 

74 

CSIRO 

G3/60 

234 

AUTUMN 

Cl 

15/11/60 

38.37 

s 

154.98 

E 

1485 

14 

75 

CSIRO 

Gl/61 

10 

SUMMER 

Cl 

18/01/61 

43.33 

S 

154.73 

E 

1799 

15 

76 

CSIRO 

G5/63 

214 

SPRING 

Cl 

13/11/63 

35.33 

s 

151.62 

E 

1500 

14 

77 

CSIRO 

G3/60 

167 

AUTUMN 

Cl 

05/08/64 

35.45 

S 

155.77 

E 

2383 

16 

78 

CSIRO 

Gl/62 

52 

SUMMER 

Cl 

14/02/62 

36.62 

S 

159.23 

E 

3406 

17 

79 

CSIRO 

Gl/61 

7 

SUMMER 

Cl 

16/01/61 

44.08 

s 

150.88 

E 

2727 

17 

80 

CSIRO 

Gl/61 

9 

SUMMER 

Cl 

17/01/61 

42.27 

s 

154.08 

E 

4019 

20 

81 

CSIRO 

Gl/61 

14 

SUMMER 

Cl 

20/01/61 

44.60 

s 

158.25 

E 

1877 

15 

82 

CSIRO 

Gl/61 

15 

SUMMER 

Cl 

20/01/61 

43.03 

s 

159.43 

E 

4146 

20 

83 

CSIRO 

Gl/61 

39 

SUMMER 

Cl 

11/01/61 

37.48 

s 

155.53 

E 

1692 

15 

84 

CSIRO 

G5/63 

197 

SPRING 

Cl 

10/11/63 

38.50 

s 

151.22 

E 

1500 

14 

85 

CSIRO 

GI/64 

74 

SUMMER 

Cl 

04/02/64 

40.  10 

s 

153.80 

E 

2500 

16 

86 

CSIRO 

Gl/61 

28 

SUMMER 

C2 

03/02/61 

35.53 

s 

166.83 

E 

1914 

16 

87 

CSIRO 

Gl/61 

25 

SUMMER 

C2 

01/02/61 

39.37 

s 

165.57 

E 

1782 

15 

88 

CSIRO 

Gl/61 

15 

SUMMER 

C2 

20/01/61 

43.03 

s 

159.43 

E 

4146 

20 

89 

CSIRO 

Gl/60 

8 

SUMMER 

C2 

04/02/60 

33.88 

s 

159.88 

E 

1948 

17 

90 

CSIRO 

Gl/60 

15 

SUMMER 

C2 

06/02/60 

35.02 

s 

167.23 

E 

2488 

16 

91 

CSIRO 

Gl/60 

17 

SUMMER 

C2 

06/02/60 

35.78 

s 

168.65 

E 

2275 

16 

92 

CSIRO 

Gl/60 

19 

SUMMER 

C2 

07/02/60 

35.88 

s 

170.42 

E 

1864 

16 

93 

CSIRO 

Gl/62 

51 

SUMMER 

C2 

13/02/62 

37.42 

s 

162.25 

E 

3431 

18 

94 

CSIRO 

Gl/62 

52 

SUMMER 

C2 

14/02/62 

36.62 

S 

159.23 

E 

3406 

17 

95 

CSIRO 

Gl/61 

16 

SUMMER 

C2 

21/01/61 

42.28 

s 

162.05 

E 

4443 

21 

96 

CSIRO 

Gl/61 

17 

SUMMER 

C2 

22/01/61 

44.85 

s 

162.80 

E 

1878 

15 

97 

CSIRO 

Gl/61 

26 

SUMMER 

C2 

02/02/61 

37.97 

s 

165.83 

E 

1594 

15 
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APPENDIX  II  CUBIC  SPLINE  COEFFICIENTS  FOR* THE  T-S  CURVES 


A  cubic  spline  fitting  procedure  fits  a  piecewise  cubic  polynomial  to  a  set  of 
data  points,  the  polynomial  coefficients  pertaining  to  adjacent  segments  being 
constrained  to  ensure  continuity  of  the  function  and  the  first  two  derivatives 
at  the  knot  points.  Two  entirely  equivalent  schemes  may  be  used  to  specify 
the  cubic  spline.  In  the  first  place,  the  polynomial  coefficients  may  be 
stated.  Alternatively,  the  knot  points  together  with  the  values  of  the 
function  and  its  first  derivative  at  the  knot  points  may  be  specified.  In 
keeping  with  the  practice  adopted  by  the  Harwell  Subroutine  Library(ref .8) ,  we 
use  the  latter  scheme. 


AREA  A1 


T  S  dS/dT 


1.7399998 

34.726685 

-0.11545753 

3 . 4487486 

34.551025 

-0.082896233 

5.1574984 

34.463531 

-0.014968872 

6.8662472 

34.503052 

0.058553696 

8.5749969 

34.642029 

0.094157219 

10.283745 

34.806015 

0.096733093 

11.992493 

34.975433 

0.10424614 

15.409996 

35.352448 

0.10685921 

18.827484 

35.620544 

0.034598351 

22.244995 

35.566254 

-0.057549477 

25.662491 

35.252518 

-0. 12745190 

29.080002 

34.667648 

-0.22149563 

AREA  A2 


T 

S 

dS/dT 

1.7799997 

34 . 699966 

-0.084634781 

5.1999969 

34.441177 

-0.023357391 

8.6199951 

34.607269 

0.096752167 

12.039993 

34.980515 

0.10944939 

15.459999 

35.330429 

0.099792480 

18.879990 

35.642120 

0.071731567 

22.299988 

35.758163 

-0.011492729 

25.719986 

35.500107 

-0.15031433 

29.139999 

34.654861 

-0.35505390 
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AREA  A3 


T  S  dS/dT 


1.8500004 

34.682465 

-0.090103149 

3.4937487 

34.517578 

-0.097715378 

5.1374979 

34.406326 

-0.022984505 

6.7812462 

34.443558 

0.054557800 

8.4249954 

34.561661 

0.088250160 

10.068741 

34.729767 

0.11484337 

11.712494 

34.925735 

0.11683750 

14.999992 

35.267654 

0.096542358 

18.287476 

35.553070 

0.069462776 

21.574982 

35.674911 

-0.0027217865 

24.862488 

35.487518 

-0.11839390 

28. 149994 

34.840240 

-0.28538132 

AREA  B1 


T  S  dS/dT 


1.1899996 

34.720779 

0.11289215 

2.7418736 

34.626480 

-0.12562466 

4.2937479 

34.494614 

-0.047614098 

5.8456221 

34.472565 

0.018503189 

7.3974972 

34.548004 

0.076808929 

10.501244 

34.888580 

0.12305641 

13.604996 

35.243362 

0. 10305786 

16.708740 

35.510086 

0.065424919 

19.812500 

35.633240 

0.012089729 

22.916245 

35.589767 

-0.036752701 

26.020004 

35.443451 

-0.048524857 

AREA  B2 


T  S  dS/dT 


1.8999996 

34.687988 

-0.072339058 

4.9149990 

34.472916 

-0.028326035 

7.9299994 

34.588760 

0.086900711 

10.944992 

34.911652 

0.11726761 

13.959999 

35.253723 

0.10568237 

16.974991 

35.530823 

0.076093674 

19.989990 

35.700470 

0.034471512 

23.004990 

35.715652 

-0.030061722 

26.020004 

35.457870 

-0.15561867 
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ARK A  B3 


T 

1.8500004 

7.8174982 

13.784996 

19.752487 

25.720001 


S 

34.691071 

34.583099 

35.247391 

35.713013 

35.447403 


AREA  Cl 


T 

1.1099997 

3.5224981 

5.9349966 

8.3474884 

10.759995 

13.172493 

15.584999 

17.997498 

20.410004 


S 

34.728088 

35.538361 

34.459381 

34.632278 

34.937805 

35.236908 

35.452240 

35.554504 

35.544281 


AREA  C2 


T  S 


1.1099997 

34.742035 

4.0724974 

34.507553 

7.0349951 

34.523254 

9.9974899 

34.843521 

12.959991 

35.205383 

15.922493 

35.453888 

18.884995 

35.580109 

21.847488 

35.633163 

24.809998 

35.646423 

dS/dT 

-0.13153839 

0.075304031 

0.11000633 

0.052697182 

-0.22023773 


dS/dT 

0.044359207 

-0.10148716 

0.027439117 

0.10853386 

0.13336277 

0.10989761 

0.066784859 

0.017894745 

-0.023929596 


dS/dT 

-0.027276993 

-0.065916061 

0.069396019 

0.12857723 

0.10705471 

0.061277390 

0.027302742 

0.011056900 

-0.0043754578 
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APPENDIX  III  DENSITY  PROFILES  FOR  THE  SOUTH-WEST  PACIFIC  OCEAN 


irt'a 

Po 

dp 

D 

n 

A1 

20.5735321 

7.50903416 

112.234772 

2.3 

A2 

20.4653168 

7.52306175 

94.4050293 

1.6 

A3 

20.3182220 

7.69532871 

98.7615967 

1.5 

B1 

23.7737579 

4.32316494 

220.479477 

2.5 

B2 

23.0042877 

5.02893066 

168.771622 

2.6 

B3 

23.6899109 

4.32462502 

194.871521 

2.5 

Cl 

25.4123688 

2.72107410 

417.001221 

2.5 

C2 

25.1103363 

2.98020172 

328.343506 

2.5 
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APPENDIX  IV  TABULATED  VALUES  OF  SOME  OCEANOGRAPHIC  VARIABLES 


IV. 1  Salinity  =  34  parts/thousand 


p 

IV. 1.1 

0°C 

Specific  volume 

5°C  10°C 

(cn^.g-1) 

15°C  20°C 

25°C 

30°C 

35°C 

40°C 

0 

.973429 

.973822 

.974496 

.975419  .976563 

.977912 

.979452 

.981172 

.983064 

100 

.968977 

.969488 

.970256 

.971252 

.972452 

.973841  .975409 

.977144 

.979040 

200 

.964662 

.965283 

.966139 

.967204 

.968457 

.969884 

.971478 

.973228 

.975129 

300 

.960477 

.961202 

.962140 

.963270 

.964572 

.966036 

.967654 

.969419 

.971325 

400 

.956417 

.957239 

.958254 

.959444 

.960793 

.962292 

.963934 

.965712 

.967623 

500 

.952476 

.953388 

.954476 

.955723 

.957116 

.958647 

.960311 

.962103 

.964018 

600 

.948651 

. 949646 

.950801 

.952101 

.953535 

.955097 

.956782 

.958586 

.960506 

700 

.944935 

.946008 

.947226 

.948576 

.950048 

.951639 

.953343 

.955158 

.957081 

800 

.941325 

.942470 

.943746 

.945142 

.946651 

.948267 

.949990 

.951816 

.953742 

900 

.937816 

.939028 

.940358 

.941797 

.943339 

.944980 

.946719 

.948554 

.950482 

1000 

.934405 

.935678 

.937057 

.938536 

.940109 

.941773 

.943527 

.945370 

.947300 

IV. 

1.2 

Isothermal 

compress ibi 

lity 

(bar- 

.1 

106) 

P 

0* 

’C 

5' 

»C 

10* 

’C 

15°C 

20‘ 

>C 

25 

°C 

30°C 

35* 

’C 

40* 

’C 

0 

46. 

46 

45 

.  18 

44. 

16 

43. 

35 

42. 

,71 

42 

.23 

41 . 

.88 

41 . 

.65 

41 . 

.52 

100 

45. 

.23 

44 

.02 

43. 

.05 

42. 

,28 

41. 

.67 

41 

.21 

40. 

.87 

40. 

.64 

40. 

.51 

200 

44. 

.05 

42 

.91 

41. 

.99 

41. 

26 

40. 

.68 

40 

.23 

39. 

.90 

39. 

,68 

39. 

.55 

300 

42 

.91 

41 

.84 

40. 

.97 

40. 

.27 

39 

.72 

39 

.29 

38. 

.97 

38 

.76 

38. 

.63 

400 

41. 

.82 

40 

.81 

39. 

.98 

39. 

32 

38, 

.79 

38 

.39 

38. 

.08 

37. 

.87 

37. 

.75 

500 

40. 

.76 

39 

.81 

39. 

.03 

38. 

.41 

37. 

.91 

37 

.52 

37, 

.23 

37, 

.03 

36. 

.91 

600 

39. 

.74 

38 

.85 

38. 

.12 

37. 

.53 

37, 

.05 

36 

.68 

36. 

.41 

36, 

.22 

36. 

10 

700 

38 

.76 

37 

.92 

37. 

.23 

36. 

,68 

36, 

.23 

35 

.88 

35. 

,62 

35. 

.44 

35. 

33 

800 

37. 

.80 

37 

.03 

36. 

.38 

35. 

.86 

35 

.43 

35 

.10 

34. 

.86 

34. 

.69 

34. 

.59 

900 

36. 

.89 

36 

.16 

35. 

.56 

35, 

,06 

34, 

.67 

34 

.36 

34, 

.  13 

33. 

.97 

33. 

.88 

1000 

36 

.00 

35 

.32 

34, 

.76 

34. 

.30 

33 

.93 

33 

.  64 

33 

.42 

33. 

.28 

33. 

.20 

J.i 
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IV. 1.3  Thermal  expansibility  (degrees- 1 . 106) 


p 

0°C 

5°C 

10°C 

15°C 

20°C 

25°C 

30°C 

35°C 

40°C 

0 

48. 

,53 

110. 

.95 

164, 

.88 

212. 

.63 

255. 

,89 

295. 

.82 

333. 

,  15 

368. 

.36 

401 

.68 

100 

76. 

12 

133. 

.  16 

182, 

.65 

226. 

.69 

266. 

78 

303. 

,95 

338. 

,86 

371 . 

.88 

403 

.22 

200 

102. 

.  14 

154. 

.  17 

199 

.54 

240. 

.10 

277. 

.21 

311. 

.80 

344. 

.42 

375, 

.39 

404 

.85 

300 

126. 

.65 

174. 

.04 

215. 

.56 

252, 

.87 

287. 

,  19 

319. 

.34 

349. 

,81 

378. 

.84 

406 

.52 

400 

149. 

,70 

192. 

.80 

230. 

.  73 

265, 

.01 

296. 

.71 

326. 

.57 

355. 

,01 

382. 

.20 

408 

.20 

500 

171. 

.34 

210. 

.48 

245 

.09 

276, 

.52 

305. 

,78 

333, 

.49 

359. 

.99 

385, 

.44 

409 

.84 

600 

191. 

.63 

227. 

.  12 

258, 

.64 

287, 

.44 

314. 

.39 

340. 

.07 

364. 

.76 

388. 

.55 

411 

.43 

700 

210. 

,62 

242. 

.75 

271 , 

.41 

297, 

.75 

322. 

55 

346. 

.32 

369. 

,28 

391. 

.50 

412 

.91 

800 

228. 

.33 

257. 

.39 

283 

.43 

307 

.47 

330. 

,26 

352. 

.22 

373. 

.56 

394. 

.27 

414 

.28 

900 

244. 

.83 

271. 

.09 

294 

.70 

316, 

,62 

337. 

,52 

357. 

.  78 

377. 

56 

396. 

,84 

415 

.50 

1000 

260. 

.  15 

283. 

.86 

305 

.25 

325, 

.20 

344. 

.33 

362, 

.99 

381. 

,30 

399 

.20 

416 

.55 

IV. 1.4  Specific  heat, 

Cp  (j -g-1 -degrees-1) 

P  0°C  5°C  10°C 

15°C  20°C  25 °C 

30°C 

35°C 

40°C 

0 

3.9926 

3.9905 

3.9918 

3.9951 

3.9991 

4.0029 

4.0059 

4.0079 

4.0090 

100 

3.9583 

3.9606 

3.9639 

3.9705 

3.9760 

3.9809 

3.9848 

3.9875 

3.9891 

200 

3.9272 

3.9333 

3.9361 

3.9479 

3.9547 

3.9605 

3.9651 

3.9683 

3.9705 

300 

3.8990 

3.9085 

3-9084 

3.9271 

3.9350 

3.9415 

3.9466 

3.9504 

3.9530 

400 

3-8735 

3.8859 

3.8809 

3.9079 

3.9167 

3.9239 

3.9294 

3.9336 

3.9366 

500 

3.8504 

3.8654 

3.8534 

3.8903 

3.8999 

3.9075 

3.9134 

3.9178 

3.9212 

600 

3.8296 

3.8468 

3.8260 

3.8742 

3.8843 

3.8923 

3.8984 

3.9031 

3.9067 

700 

3.8108 

3.8299 

3.7986 

3.8594 

3.8700 

3.8782 

3.8845 

3.8893 

3.8932 

800 

3.7939 

3.8146 

3.7714 

3.8458 

3.8567 

3.8651 

3.8715 

3.8765 

3.8806 

900 

3.7786 

3.8008 

3.7443 

3.8334 

3.8445 

3.8530 

3.8594 

3.8645 

3.8687 

000 

3.7649 

3.7883 

3.7173 

3.8221 

3.8334 

3.8418 

3.8482 

3.8534 

3.8577 
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IV. 1.5  Sound  speed  (m.s-1) 


P  0°C  5°C  10°C  15°C  20°C  25°C  30°C  35°C  40 °C 


0  1447.80  1447.80  1447.80  1447.80  1447.80  1447.80  1447.80  1447.80  1447.80 

100  1464.19  1464.19  1464.19  1464.19  1464.19  1464.19  1464.19  1464.19  1464.19 

200  1480.90  1480.90  1480.90  1480.90  1480.90  1480.90  1480.90  1480.90  1480.90 

300  1497.90  1497.90  1497.90  1497.90  1497.90  1497.90  1497.90  1497.90  1497.90 

400  1515.15  1515.15  1515.15  1515.15  1515.15  1515.15  1515.15  1515.15  1515.15 

500  1532.62  1532.62  1532.62  1532.62  1532.62  1532.62  1532.62  1532.62  1532.62 

600  1550.26  1550.26  1550.26  1550.26  1550.26  1550.26  1550.26  1550.26  1550.26 

700  1568.05  1568.05  1568.05  1568.05  1568.05  1568.05  1568.05  1568.05  1568.05 

800  1585.94  1585.94  1585.94  1585.94  1585.94  1585.94  1585.94  1585.94  1585.94 

900  1603.90  1603.90  1603.90  1603.90  1603.90  1603.90  1603.90  1603.90  1603.90 

1000  1621.90  1621.90  1621.90  1621.90  1621.90  1621.90  1621.90  1621.90  1621.90 


IV. 2  SaLinity  =  35  parts/thousand 


IV. 2.1  Specific  volume  (cn^.g-1) 


P  0°C  5°C  10°C  15°C  20°C  25°C  30°C  35°C  40°C 


0  .972664  .973070  .973755  .974685  .975835  .977189  .978733  .980456  .982350 
100  .968227  .968750  .969527  .970530  .971735  .973130  .974700  .976438  .978336 
200  .963926  .964558  .965423  .966494  .967751  .969183  .970779  .972532  .974435 
300  .959755  .960489  .961436  .962570  .963877  .965345  .966965  .968732  .970640 
400  .955707  .956538  .957560  .958755  .960108  .961610  .963254  .965035  .966947 
500  .951779  .952699  .953793  .955044  .956440  .957974  .959640  .961434  .963350 
600  .947965  .948968  .950128  .951433  .952869  .954434  .956120  .957925  .959846 
700  .944261  .945340  .946563  .947916  .949391  .950984  .952689  .954506  .956430 
800  .940661  .941812  .943093  .944492  .946002  .947620  .949344  .951170  .953097 
900  .937163  .938379  .939713  .941155  .942698  .944341  .946080  .947916  .949845 
1000  .933761  .935038  .936421  .937902  .939476  .941141  .942896  .944739  .946670 
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IV. 2. 2 

Isothermal  compressibility  (bar- 1 . 10e) 

p 

0°C 

5°C  10°C  15°C  20°C  25°C 

30°i: 

35°C 

40° 

0 

46 . 

.34 

45. 

.07 

44. 

.06 

43. 

.25 

42 

.62 

42. 

14 

41 

.80 

41 

.57 

41 . 

.45 

100 

45. 

.  1  1 

43. 

.92 

42. 

.96 

42. 

.  19 

41  . 

.59 

41 . 

13 

40. 

.  79 

40. 

.57 

40. 

44 

200 

43. 

.94 

42. 

.81 

41  . 

.90 

41 . 

.  17 

40. 

.59 

40. 

15 

39. 

.83 

39. 

.61 

39. 

.48 

300 

42. 

.81 

41 . 

.74 

40. 

.88 

40 

.  19 

39, 

.  64 

39. 

22 

38. 

.90 

38, 

.69 

38. 

.56 

400 

41 . 

,72 

40 

.71 

39. 

.90 

39 

.24 

38. 

.  72 

38. 

32 

38. 

.02 

37. 

.81 

37. 

.69 

500 

40. 

,66 

39. 

.  72 

38. 

.95 

38. 

.33 

37. 

.83 

37. 

45 

37. 

.  16 

36. 

.96 

36. 

.85 

600 

39. 

.65 

38. 

.76 

38 

.04 

37 

.45 

36, 

.98 

36. 

62 

36. 

.34 

36. 

.  16 

36. 

.04 

700 

38. 

.67 

37, 

.84 

37 

.  16 

36. 

.61 

36 

.16 

35. 

82 

35, 

.56 

35, 

.38 

35. 

,27 

800 

37. 

.72 

36 

.95 

36. 

.31 

35 

.79 

35. 

.37 

35. 

04 

34. 

.80 

34. 

.63 

34. 

.54 

900 

36. 

.81 

36, 

.09 

35, 

.49 

35, 

.00 

34. 

.61 

34. 

30 

34. 

.07 

33. 

.92 

33. 

.83 

1000 

35. 

.93 

35 

.25 

34 

.70 

34 

.24 

33, 

.87 

33. 

58 

33, 

.37 

33. 

.23 

33. 

.  15 

IV, 

,2.3 

Thermal 

expans i 

[bill 

ity 

(degrees- 

-1 .  10 

i6) 

P 

O' 

aC 

5°C 

10°C 

15°C 

2  0°C 

25  °C 

30  °C 

35°C 

40°C 

0 

51 

.71 

113. 

.46 

166. 

.85 

214. 

.  18 

257. 

.13 

296. 

.84 

334. 

.05 

369. 

.20 

402. 

.54 

100 

79 

.05 

135. 

.47 

184, 

.47 

228. 

,  12 

267. 

92 

304. 

.89 

339. 

,69 

372. 

.68 

404. 

.04 

200 

104 

.84 

156. 

,30 

201 . 

.21 

241  . 

.41 

278. 

.26 

312. 

.67 

345. 

.  19 

376. 

.  13 

405. 

.63 

300 

129 

.  12 

176. 

.00 

217. 

.09 

254. 

.07 

288. 

15 

320. 

.  14 

350. 

.52 

379. 

.53 

407. 

.27 

400 

151 

.97 

194. 

59 

232. 

.  13 

266. 

,  10 

297. 

58 

327. 

,30 

355. 

.66 

382. 

.85 

408. 

90 

500 

173 

.42 

212. 

,  12 

246. 

,36 

277. 

,52 

306. 

57 

334. 

,  15 

360. 

.60 

386. 

.05 

410. 

.51 

600 

193 

.54 

228. 

.61 

259. 

.80 

288. 

.33 

315. 

,  10 

340. 

.  66 

365. 

.30 

389. 

.11 

412. 

.05 

700 

212 

.36 

244. 

.  11 

272. 

.47 

298. 

.56 

323. 

19 

346. 

.85 

369. 

.78 

392. 

,01 

413. 

,50 

800 

229 

.93 

258. 

.  64 

284. 

.38 

308. 

,20 

330. 

82 

352. 

.69 

373. 

.99 

394. 

,74 

414. 

82 

900 

246 

.30 

272 

.22 

295 

.56 

317, 

.27 

338. 

02 

358. 

.  19 

377. 

.95 

397. 

.26 

415. 

.99 

000 

261 

.51 

284 

.90 

306. 

.03 

325. 

.78 

344. 

77 

363. 

.35 

381 . 

.63 

399, 

.57 

417. 

.00 
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IV. 2. 4 

Specific  heat 

p 

0°C 

5°C 

10°C 

0 

3.9865 

3.9847 

3.9864 

100 

3.9526 

3.9552 

3.9587 

200 

3.9219 

3.9282 

3.9312 

300 

3.8940 

3.9036 

3.9038 

400 

3.8688 

3.8813 

3.8765 

500 

3 . 8400 

3.8010 

3.8493 

000 

3.8254 

3.8420 

3.8222 

700 

3.8069 

3.8200 

3.7951 

800 

3.7902 

3.8109 

3.7682 

900 

3.7751 

3.7973 

3.7413 

1000 

3.7615 

3.7849 

3.7146 

IV. 2. 5 

Sound 

speed  (i 

P 

0°C 

5°C 

10°C 

,  Cp  (j .g-1 .degrees-1) 


15°C 

20°C 

25°C 

3.9899 

3.9940 

3.9980 

3.9655 

3.9712 

3.9762 

3.9431 

3.9500 

3.9559 

3.9224 

3.9304 

3.9370 

3.9035 

3.9123 

3.9194 

3.8860 

3.8955 

3.9031 

3.8700 

3.8801 

3.8879 

3.8553 

3.8658 

3.8738 

3.8419 

3.8526 

3.8608 

3.8296 

3.8405 

3.8488 

3.8184 

3.8294 

3.8376 

.s-1) 

15°C  20°C  25°C 


30°C 

35°C 

40QC 

4.0011 

4.0032 

4.0043 

3.9801 

3.9828 

3.9844 

3.9604 

3.9637 

3.9658 

3.9420 

3.9457 

3.9482 

3.9249 

3.9289 

3.9318 

3.9089 

3.9132 

3.9163 

3.8939 

3.8984 

3.9019 

3.8800 

3.8846 

3.8883 

3.8670 

3.8718 

3.8756 

3.8549 

3.8598 

3.8637 

3.8437 

3.8486 

3.8526 

30°C  35°C  40°C 


0 

100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 


1449. 14  1449.14  1449.14 

1465.54  1465.54  1465.54 

1482.26  1482.26  1482.26 

1499.27  1499.27  1499.27 

1516.51  1516.51  1516.51 

1533.97  1533.97  1533.97 

1551.60  1551.60  1551.60 

1569.37  1569.37  1569.37 

1587.24  1587.24  1587.24 

1605.18  1605.18  1605.18 

1623.15  1623.15  1623.15 


1449.14  1449.14  1449.14 

1465.54  1465.54  1465.54 

1482.26  1482.26  1482.26 

1499.27  1499.27  1499.27 

1516.51  1516.51  1516.51 

1533.97  1533.97  1533.97 

1551.60  1551.60  1551.60 

1569.37  1569.37  1569.37 

1587.24  1587.24  1587.24 

1605.18  1605.18  1605.18 

1623.15  1623.15  1623.15 


1449.14  1449.14  1449.14 

1465.54  1465.54  1465.54 

1482.26  1482.26  1482.26 

1499.27  1499.27  1499.27 

1516.51  1516.51  1516.51 

1533.97  1533.97  1533.97 

1551.60  1551.60  1551.60 

1569.37  1569.37  1569.37 

1587.24  1587.24  1587.24 

1605.18  1605.18  1605.18 

1623.15  1623.15  1623.15 


Li 


J 
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IV.  3  S.tlitiily  =  36  parts/thousand 
IV.  5.1  Specific  volume  (cm3.g-') 


P  0oC  5°C  10°C  15°C  20°C  25°C  30°C  35°C  4o°C 


0  .971900  .972320  .973015  .973953  .975109  .976467  .978013  .979739  .981637 

100  .967478  .968014  .968800  .969810  .971021  .972418  .973991  .975732  .977634 

200  .963190  .963835  .964708  .965785  .967047  .968482  .970080  .971836  .973742 

300  .959033  .959779  .960732  .961872  .963183  .964654  .966276  .968046  .969956 

400  .954998  .955839  .956868  .958068  .959425  .960929  .962574  .964357  .966272 

500  .951082  .952011  .953111  .954367  .955766  .957302  .958969  .960764  .962684 

600  .947280  .948291  .949457  .950765  .952205  .953770  .955457  .957264  .959187 

700  .943587  .944674  .945901  .947258  .948735  .950328  .952035  .953852  .955779 

800  .939998  .941156  .942440  .943842  .945354  .946973  .948697  .950525  .952454 

900  .936509  .937732  .939069  .940513  .942059  .943702  .945441  .947278  .949209 

1000  .933117  .934400  .935786  .937269  .938844  .940509  .942264  .944108  .946040 


IV. 3. 2  Isothermal  compressibility  (bar_i.106) 


p 

0‘ 

5C 

5°C 

10°C 

15°C 

20°C 

25°C 

30°C 

35‘ 

JC 

40°C 

0 

46. 

.22 

44. 

.96 

43. 

96 

43. 

16 

42. 

.54 

42. 

06 

41. 

,72 

41 . 

49 

41 . 

.37 

100 

45. 

.00 

43. 

.82 

42. 

.86 

42. 

,  10 

41 . 

.50 

41 . 

.05 

40. 

.72 

40. 

,49 

40. 

.37 

200 

43. 

.83 

42. 

.71 

41 . 

,81 

41 . 

08 

40. 

.51 

40. 

08 

39. 

.75 

39. 

54 

39. 

.41 

300 

42. 

.70 

41 . 

.65 

40. 

.79 

40. 

,  11 

39. 

.56 

39. 

.  14 

38 

.83 

38. 

.62 

38. 

.50 

400 

41 . 

.61 

40. 

.62 

39. 

.81 

39. 

.  16 

38, 

.65 

38. 

.25 

37. 

.95 

37. 

,74 

37. 

.62 

500 

40 

.57 

39. 

.63 

38 

.87 

38. 

.25 

37 

.76 

37. 

.38 

37. 

.  10 

36. 

.90 

36. 

.79 

600 

39. 

.56 

38. 

.68 

37 

.96 

37. 

,38 

36. 

.91 

36. 

.55 

36 

.28 

36. 

.  10 

35, 

.99 

700 

38 

.58 

37. 

.76 

37 

.08 

36, 

.54 

36 

.  10 

35. 

.75 

35 

.50 

35. 

.32 

35 

.22 

800 

37 

.64 

36. 

.87 

36 

.24 

35. 

,  72 

35. 

.33 

34. 

.98 

34. 

.74 

34. 

.58 

34, 

.48 

900 

36 

.73 

36. 

.01 

35. 

.42 

34. 

,94 

34. 

.55 

34. 

.24 

34. 

.02 

33. 

.86 

33 

.78 

1000 

35 

.85 

35 

.  19 

34. 

.63 

34. 

.  18 

33, 

.81 

33. 

.53 

33 

.32 

33. 

.  18 

33 

.  10 
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IV. 3. 3  Thermal  expansibility  (degrees-1 . 106) 


p 

0°C 

5°C 

10°C 

15°C 

20°C 

25°C 

30°C 

35°C 

40°C 

0 

54. 

.89 

115, 

.96 

168 

.80 

215. 

.71 

258, 

.35 

297. 

.85 

334. 

,94 

370. 

.06 

403. 

.43 

100 

81. 

.98 

137, 

.78 

186 

.27 

229. 

.53 

269 

.04 

305, 

.83 

340. 

.53 

373, 

.49 

404. 

.90 

200 

107. 

.54 

158. 

.43 

202 

.86 

242. 

.70 

279, 

.29 

313. 

.53 

345. 

,97 

376. 

.90 

406. 

.45 

300 

131 

.60 

177. 

.95 

218 

.61 

255. 

.25 

289. 

.09 

320, 

.93 

351. 

.24 

380. 

.25 

408. 

.05 

400 

154. 

.24 

196. 

.38 

233, 

.52 

267. 

.  18 

298, 

.44 

328. 

.02 

356. 

.32 

383. 

.52 

409. 

.64 

500 

175, 

.51 

213. 

.  75 

247 

.63 

278. 

.50 

307 

.35 

334, 

.80 

361 . 

,20 

386. 

.67 

411 , 

.21 

600 

195. 

.45 

230 

.  1 1 

260 

.95 

289, 

.22 

315 

.80 

341 , 

.25 

365. 

.85 

389. 

.69 

412, 

.71 

700 

214. 

.  1 1 

245. 

.47 

273, 

.51 

299. 

.35 

323. 

.81 

347, 

.37 

370. 

.27 

392. 

.54 

414. 

.  1 1 

800 

231 

.54 

259. 

.88 

285 

.32 

308, 

.91 

331 

.38 

353. 

.  16 

374, 

.44 

395. 

.22 

415 

.39 

900 

247. 

.  78 

273. 

.36 

296. 

.41 

317. 

.91 

338. 

.51 

358. 

,60 

378. 

.34 

397. 

,69 

416. 

.52 

1000 

262, 

.87 

285. 

.94 

306 

.80 

326. 

.34 

345, 

.20 

363. 

.70 

381. 

.97 

399. 

.95 

417. 

.48 

IV. 3. 4  Specific  heat, 

Cp  (j -g-1. degrees-1) 

P  0°C  5°C  10°C 

15°C  20°C  25°C 

30°C 

35°C 

40°C 

0 

3.9804 

3.9790 

3.9809 

3.9847 

3.9890 

3.9931 

3.9964 

3.9985 

3.9996 

100 

3.9470 

3.9498 

3.9536 

3.9605 

3.9663 

3.9714 

3.9754 

3.9781 

3.9797 

200 

3.9166 

3.9231 

3.9264 

3.9383 

3.9453 

3.9512 

3.9558 

3.9590 

3.9610 

300 

3 . 889 1 

3.8988 

3.8992 

3.9178 

3.9258 

3.9324 

3.9374 

3.9410 

3.9435 

400 

3.8642 

3.8768 

3.8722 

3.8990 

3.9078 

3.9149 

3.9203 

3.9242 

3.9270 

500 

3.8417 

3.8567 

3.8453 

3.8817 

3.8912 

3.8986 

3.9043 

3.9084 

3.9115 

600 

3.8214 

3.8386 

3.8184 

3.8658 

3.8758 

3.8835 

3.8894 

3.8937 

3.8969 

700 

3.8030 

3.8221 

3.7916 

3.8513 

3.8616 

3.8695 

3.8754 

3.8799 

3.8833 

800 

3.7865 

3.8072 

3.7650 

3.8380 

3.8485 

3.8565 

3.8625 

3.8670 

3.8705 

900 

3.7717 

3.7937 

3.7384 

3.8258 

3.8365 

3.8445 

3.8504 

3.8549 

3.8586 

000 

3.7583 

3.7816 

3.7119 

3.8147 

3.8255 

3.8334 

3.8392 

3.8437 

3.8474 
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IV. 3. 5  Sound  speed 


P  0°C  5°C  10 


0  1450.48  1450.48  1450 
100  1466.90  1466.90  1466 
200  1483.63  1483.63  1483 
300  1500.63  1500.63  1500 
400  1517.87  1517.87  1517 
500  1535.32  1535.32  1535 
600  1552.93  1552.93  1552 
700  1570.68  1570.68  1570 
800  1588.53  1588.53  1588 
900  1606.45  1606.45  1606 
1000  1624.40  1624.40  1624 


(m . s- 1 ) 


C  15°C  20°C  25 


48  1450.48  1450.48  1450 
90  1466.90  1466.90  1466 
63  1483.63  1483.63  1483 
63  1500.63  1500.63  1500 
87  1517.87  1517.87  1517 
32  1535.32  1535.32  1535 
93  1552.93  1552.93  1552 
68  1570.68  1570.68  1570 
53  1588.53  1588.53  1588 
45  1606.45  1606.45  1606 
40  1624.40  1624.40  1624 


C  30°C  35 °C  40°C 


48  1450.48  1450.48  1450.48 
90  1466.90  1466.90  1466.90 
63  1483.63  1483.63  1483.63 
63  1500.63  1500.63  1500.63 
87  1517.87  1517.87  1517.87 
32  1535.32  1535.32  1535.32 
93  1552.93  1552.93  1552.93 
68  1570.68  1570.68  1570.68 
53  1588.53  1588.53  1588.53 
45  1606.45  1606.45  1606.45 
40  1624.40  1624.40  1624.40 


o  o  n  o 
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APPENDIX  V  LISTING  OF  A  FORTRAN  IV  PROGRAMME  TO  PROCESS  XBT  DATA 


V. 1  Main  programme 


C  XBT. PROFILE. FORT.  A  PROGRAMME  TO  CARRY  OUT  PRE-PROCESSING  OF 
C  XBT  DATA. 

C 

C  THE  STANDARD  OUTPUT  FROM  THIS  PROGRAMME  IS  A  TABLE  OF 
C  SIGNIFICANT  OCEANOGRAPHIC  PARAMETERS. 

C  PLOTTING  OPTIONS  ARE  SELECTED  FROM  STANDARD  INPUT  DEVICE 
C  XBT  DATA  IS  ASSUMED  AVAILABLE  FROM  LOGICAL  UNIT  2 

C  IF  LNGPLT  IS  SELECTED  A  DATASET  IS  ASSUMED  AVAILABLE  ON  LOGICAL  UNIT  3 
C  FOR  WRITING  DATA  IN  PREPARATION  FOR  PLOTTING. 

REAL* 8  KO(110),A(110),B(110),DER(110),W(330),X8(110),Y8(110) 

REAL*4  CRUISE (9) ,SDEP(2) ,BPRESS(2) ,SV(2) ,UPI (2) ,LAT, 

1  LONG , VP ( 1 10) , TITLE 1 (5) ,TITLE2(5) 

REAL*4  Z(110),T(110),S(110) ,DHA( 110) ,DELTA( 1 10) ,RHO( 1 10) ,KP( 1 10) , 

1  SIGMAT(llO) ,ALPHA(110) , BETA(llO) ,V(110),P(110),G(110),V0(110), 

2  CP(110) 

DATA  TITLE  1/' TEMP' ,'ERAT' ,'URE  DE’.'G.  C’/ 

DATA  TITLE2/ ' SOUN' , ' D  SP' , 

1  '  EED  VM/S  V  '  /  ,SV(  1 )  ,SV(2)/0.  , 1480 . / ,UPI ( 1 ) ,UPI (2)/ 10.  , 

2  20./ 

******* 

CHOOSE  PLOT  OPTIONS 


READ(5 ,*)  IPLT, ISCALE , ILPLT, ITYP 
WRITE (6, 350)  IPLT, ISCALE, ILPLT, ITYP 

C 

C  IPLT  =  1  PLOT  SQUARE  PLOT 
C  ISCALE  =  1  FOR  SHALLOW  SQUARE  PLOT 
C  =  2  FOR  DEEP  SQUARE  PLOT 

C  ILPLT  =  1  WRITES  OUTPUT  TO  DISC  READY  FOR  LONG  PLOT 
C  ITYP  =  1  TEMPERATURE  LONG  PLOT 
C  =2  VELOCITY  LONG  PLOT 

C  =3  BOTH  LONG  PLOTS 

C  READ  HOW  MANY  XBTS  IN  BATCH 
READ(2 , 10)  MANY 
10  FORMAT(IIO) 

C  LOOP  ON  ENTIRE  PROGRAM  'MANY'  TIMES 
DO  90  MM=1 .MANY 
READ (2 , 190) (CRUISE(I) ,1=1,9) 

READ( 2 , 200 ) DY , XMTH , YR , TIME , LONG , LAT 

READ ( 2 , 2 40 ) SDEP ( 1 ) , SDEP ( 2 ) , WETT , DRYT , WD I R , WSP , BPRESS ( 1 ) , 

Z  BPRESS (2) 

READ(2,250)IJ,L 

IF( ILPLT  .EQ.l)  WRITE(3,250)IJ,L 
C  EVALUATE  SURFACE  VALUE  OF  GRAVITATIONAL  ACCELERATION 
RLAT  =  0.01 745 32925*LAT 

GO  =  9. 78049*(1 .0  +  0.0052884*(SIN(RLAT))**2  -  0.0000059* 
1  (SIN(2 . 0*RLAT) )**2) 


m  J 


4J 
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C  /.  STOKED  POSITIVE  DOWNWARDS 
i:  RKA1)  IN  TEMPERATURE  PROFILE. 

REAL)  (2 ,  320) (Z( I ) ,T( I ) , 1=1  ,  L) 

C  COMPUTE  SALINITY  AND  SIGMA-T  PROFILES. 

C  COMMENCE  LOOP  ON  I. 

00  20  1=1, L 

C  EVALUATE  LOCAL  VALUE  OF  GRAVITATIONAL  ACCELERATION. 

G(I )  =  GO* (1.0  +  2 . 26E-07*Z ( I ) ) 

C  COMPUTE  SALINITY  ( PARTS/ THOUSAND ) . 

S(I)  =  SAL ( LAT , LONG ,  T  ( I )  ) 

CALL  DSWI(SCr),T(r) , VO(I) ,RHO(I) , SIGMAT(I)) 

CALL  DSW1A(S(I) ,T(I),KO(I),A(I),B(I)) 

C  AS  AN  INITIAL  APPROXIMATION,  SET  SPECIFIC  VOLUME  UNDER  PRESSURE  EQUAL 
C  TO  VALUE  AT  ATMOSPHERIC  PRESSURE. 

VP(I)  =  VO(I) 

C  CONVERT  Z  AND  RHO  TO  DOUBLE  PRECISION 
Y8(I)  =  DBL£(Z(I)) 

X8(I )  =  DBLE(RHO( I ) ) 

C  END  OF  LOOP  ON  I. 

20  CONTINUE 

C  COMMENCE  ITERATION  TO  CORRECT  FOR  PRESSURE  BY  FITTING  CUBIC  SPLINE 
C  TO  CURRENT  DENS  TTY  ESTIMATES. 

PCD  =  o.o 
30  CONTINUE 

CALL  TB04AD(L,Y8,X8,DER,W) 

IF  (W(l).EQ.l.O)  WRITE(6,40) 

40  FORMAT ('  UNSUCCESSFUL  RETURN  FROM  TB04A. 1 ) 

DV  =  0.0 

DO  50  J=2,L 

C  LOOP  ON  DATA  POINTS  TO  COMPUTE  PRESSURES  AND  SUBSEQUENTLY  DENSITY 
C  CORRECTIONS  (INDEX  J) . 

P(J)  =  1.0E-02*G(J)*QG02AD(0.0,Y8(J) ,L,Y8,X8,DER) 

CALL  DSW2(P(J) ,V0(J) ,K0(J) ,A(J) ,B(J) ,KP(J) ,VN,BETA(J) ) 

IF  (  ABS (VN-VP ( J ) ) . GT . DV)  DV  =  ABS(VN  -  VP(J)) 

VP(J)  =  VN 
RHO(J)  =  1.0/VP(J) 

X8(I)  =  DBLE(RHO(J) ) 

C  END  OF  LOOP  ON  J. 

50  CONTINUE 

C  CHECK  WHETHER  FURTHER  ITERATIONS  ARE  REQUIRED. 

IF  (DV.GT.1.0E-08)  GO  TO  30 

C  EVALUATE  SPECIFIC  VOLUME  ANOMALY. 

C  SAL  AT  35PPT  TEMP  AT  ODEG 

CALL  DSW1(35. ,0. ,VOF,RHOF,STF) 

CALL  DSW1A(35. ,0. ,KOF,AF,BF) 

DO  60  1=1, L 

CALL  DSW2 ( P ( I ) , VOF , KOF , AF , BF , KPF , V2 , BETAF ) 

C  COMPUTE  SPECIFIC  VOLUME  ANOMALY. 

C  END  OF  LOOP  ON  I. 

DELTA(I)  =VP(I)-  V2 

C  COMPUTE  SOUND  SPEED  FROM  TEMP, SAL,  AND  PRESSURE. 

CALL  SSSW1(P(I),S(I),T(I),V(I)) 

C  COMPUTE  EXPANSIBILITY  AND  SPECIFIC  HEAT. 

CALL  PSW1 (S(I),T(I),P(I),V0(I),VP(I),K0(I),A(I),B(1) , ALPHA ( I 
1  ),  CP(I)) 

C  CONVERT  P  AND  DELTA  TO  DOUBLE  PRECISION 
Y8( I )  =  DBLE(P( I ) ) 
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X8(  I )  =  DBLE (DELTA ( l ) ) 

60  CONTINUE 

C  FIT  A  CUBIC  SPLINE  TO  SPECIFIC  VOLUME  ANOMALY  VS.  PRESSURE. 

CALL  TB04AIX L,Y8,X8,DER,W) 

IF  (W(  I  )  .EQ.  I  .0)  WR 1 TF, ( (> , 40 ) 

C  EVALUATE  DYNAMIC  HEIGHT  ANOMALY,  REFERRED  TO  THE  SURFACE. 

DO  70  1=1, L 

70  DHA(I)  =  10.*QG02AD(0.0,Y8(I),L,Y8,X8,DER) 

C  OUTPUT  PROFILE. 

WRITE (6, 220) (CRUISE(I) ,1=1,9) 

WRITE(6 , 160)  IJ 
WRITE(6 , 170) 

WRITE(6 , 180) 

WRITE(6 ,210)  DY,XMTH,YR, TIME, LONG, LAT 
WRITE(6 , 170) 

WRITE (6 , 230) 

WRITE(6 , 260)SDEP( 1 ) ,SDEP(2) ,WETT ,DRYT , WDIR , WSP ,BPRESS( 1 ) , 

Z  BPRESS(2) 

WRITE(6 , 170) 

WRITE(6 , 270) 

WRITE  16,280) 

WRITE(6 ,290) 

WRITE(6,300) 

WRITE(6 ,310) 

WRITE(6 ,270) 

WR I TE ( 6 , 1 7  0 ) 

DO  80  1  =  1,  L 

ALPHA II)  =  ALPHA(I)  *  1E06 
BETA ( 1 )  =  BETA(I)  *  1E06 

WRITE  lb, 330)  Z(I),T(I),S(I) ,DHA( I ) ,DELTA(I ) ,RHO(I ) , 

Z  SIGMAT( I ) , ALPHA ( I ) ,BETA(I ) ,V(I),P(I),CP(I) 

IF( ILPLT.EQ. 1) 

Y  WRITE (3, 3 AO)  Z(I ) ,T(I ) ,S( I ) ,DHA(I ) ,DELTA(I ) ,RHO(I ) , 

Z  SIGMAT( I ) ,ALPHA(I ) ,BETA(I ) ,V(I),P(I),CP(I) 

80  CONTINUE 

WRITE (6, 170) 

IF(IPLT.EQ.l)  CALL  PLOT1 (Z,T,S,SIGMAT,V,L, IJ.ISCALE) 

>)0  CONTINUE 

lF(lLPLT.NE.l)  GO  TO  150 

rF(ITYP.EQ. 3)  NPLT  =  2 
TF( ITYP .NE . 3)  NPLT  =  1 
DO  140  I N= 1 ,NPLT 
REWIND  3 

DO  130  IK=l .MANY 
READ (3, 250)  IJ,L 
DO  100  1=1, L 

READ ( 3 , 360 )  Z( I ) ,T( I ) ,V( I ) 

100  CONTINUE 

IF(lTYP.EQ.l)  GO  TO  110 
IF(ITYP.EQ.2)  GO  TO  120 
IF(ITYP.EQ.3.AND.IN.EQ.l)  GO  TO  110 
IF ( ITYP . EQ. 3 .AND. IN .EQ. 2)  GO  TO  120 
110  IF(IK.EQ.l)  CALL  LNGPLT(SV(1) ,UPI(1) ,TITLE1 ,T,Z,L, 

1  IJ, IK, MANY) 

IF(IK.NE.l)  CALL  ADDPLT(T,Z,L,IJ,IK) 

GO  TO  130 

120  IF(IK.F.Q.l)  CALL  LNGPLT(SV(2)  ,UPI  (2)  ,TITLE2,V,Z,L, 
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1  IJ, IK, MANY) 

IF(IK.NE.l)  CALL  ADDPLT(V,Z,L,IJ,IK) 

130  CONTINUE 

140  CONTINUE 
CALL  ENDGR 
150  STOP 

160  FORMAT (47X, 'XBT  NUMBER', 13) 

170  FORMAT ( 1 HO) 

180  FORMAT (32X, 'DATE' ,7X,'TIME' ,6X, ' LONGTITUDE' ,8X, 'LATITUDE' ) 

190  FORMAT ( 9 A4) 

200  FORMAT (A2 , IX, A2 , IX , A2 , 3X ,A4,4X,F7 . 2 , 1 1X,F6 . 2) 

210  FORMAT (30X,2(A2, '/’ ) ,A2,5X,A4,5X,F6.2, '  DEG.E ’ ,6X,F5 . 2 , ’  DEG. S’) 
220  F0RMAT(1H1 ,37X,9A4) 

230  FORMAT (17X, 'SONIC  DEPTH' , 7X, 'AIR  TEMPERATURE (DEG. C) ', 1 IX, ' WIND' , 
Z12X, 'ATMOS. PRESS ' ) 

240  FORMAT ( 1 X , A3 , A2 , 1 OX , A4 , 5X , A4 , 6X , A3 , A3 , 9X , A4 , A2 ) 

250  FORMAT(I3,I5) 

260  FORMAT (19X, A3, A2, 'M' , 10X, 'WET  ',A4,5X,'DRY  ',A4,5X,’DIR  ’, 
ZA3,3X,A3, '  KNOTS' ,5X,A4,A2,'  MB') 

270  FORMAT ( ' - - - - 

. . . . . — - - ') 

280  FORMAT ( '  DEPTH  TEMPERATURE  SALINITY  DYNAMIC  HT.  DELTA  DENSITY 
Z  SIGMA-T  ALPHA  BETA  SOUND  PRESSURE  SPECIFIC’) 

290  FORMAT (3 IX, 'ANOMALY' ,53X'SPEED' ,16X,'HEAT’) 

300  FORMAT (2X, 'M' ,9X, 'C' ,8X, 'PPT' ,7X, 'DYNM' ,6X'CM**3/GM' ,2X 

Z, 'GM/CM**3' , '  (GM/CM**3) ' ,2X, '*10**6/' ,2X, '*10** 6/ ' ,3X, 'M/SEC' ,6X, 
Y  '  BAR') 

310  FORMAT (62X, '/ 1000' ,6X,'DEG.C' ,5X,’BAR' ,24X, ' J/G/DEG' ) 

320  FORMAT (12F6. 2) 

330  FORMAT ( IX, F6 .0,4X,F5 . 2 ,5X,F5 . 2 ,5X ,F6 . 3 , 4X,F8 . 6 ,2X,F8 . 6 , 2X,F7 . 3 , 
Z5X,F5.1,4X,F5.2,3X,F7.2,4X,F6.3,4X,F6.4) 

340  FORMAT (F6 . 0 ,F5 . 2 ,F5 . 2 ,F6 . 3 ,F8 . 6 ,F8 . 6 ,F7 . 3 , 

ZF5 . 1 ,F5 . 2 ,F7 . 2 ,F7 . 1 ,F6 . 4) 

350  FORMAT ( '  IPLT  =  ’,  12,’  ISCALE  =  ’,12,’  ILPLT  =  ',  12,  ’ ITYP  =  ’, 

Z  12) 

360  FORMAT (F6 . 0 , F5 . 2 , 44X , F7 . 2 ) 

END 


V.2  Subroutine  DSW1 


SUBROUTINE  DSW1(S,T,V0,R0,SIGMAT) 

C  A  SUBROUTINE  TO  DETERMINE  THE  SPECIFIC  VOLUME  OF  SEAWATER  AT 
C  ATMOSPHERIC  PRESSURE  AND  SPECIFIED  SALINITY  AND  TEMPERATURE  USING  THE 
C  SOUND-SPEED  DERIVED  EQUATION  OF  STATE  OF  CHEN  AND  MIlLERO  (1). 

C  INPUT  PARAMETERS: 

C  S  REAL*4 .  SALINITY  (PARTS/THOUSAND). 

C  T  REAL*4 .  TEMPERATURE  (DEGREES  CENTIGRADE). 

C  OUTPUT  PARAMETERS: 

C  VO  REAL*4 .  SPECIFIC  VOLUME  OF  SEAWATER  (CM**3/G)  AT  ATMOSPHERIC 
C  PRESSURE  AND  SPECIFIED  SALINITY  AND  TEMPERATURE. 

C  RO  REAL*4 .  CORRESPONDING  DENSITY  (G/CM**3) . 

C  SIGMAT  REAL*4 .  CORRESPONDING  VALUE  OF  SIGMA-T  ( (G/CM**3)*10**-3) . 

C  REFERENCE: 


til 
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C  (1)  CHEN,  C.-T.,  AND  MILLERO,  F.J.  1978.  THE  EQUATION  OF  STATE  OF 
C  SEAWATER  DETERMINED  FROM  SOUND  SPEEDS.  J.  MAR.  RES.,  VOL.  36, 

C  657-691. 

C 

REAL*8  AO , A 1 , A2 , A3 , A4 , A5 , A6 
REAL*8  BO.Bl ,B2,B3,B4 
REAL*8  CO , Cl ,C2 ,C3 
REAL*8  DO,Dl ,D2,D3 
REAL* 8  AT,BT,CT,DT 

DATA  AO , A 1 , A2 , A3 , A4 , A5 , A6/0 . 9998395 , 6 . 79 14D-05 , -9 . 0894D-06 , 

1  1 .0171D-07.-1 .28460-09,1 . 1592D-1 1 , -5 . 0125D- 14/ 

DATA  BO , B 1 , B2 , B3 , B4/8 . 259 1 7D-04 , -4 . 4490D-06 , 1 . 0485D-07 , 

1  -1.25800-09,3.3150-12/ 

DATA  CO, Cl , C2 , C3/-6 . 33761D-06 ,2 . 8441D-07 , -1 . 687 ID-08 ,2 . 83258D-10/ 
DATA  DO ,D1 ,D2 ,D3/5 . 4705D-07 , -1 . 97975D-08 , 1 . 6641D-09 , -3 . 1203D-1 1/ 

C  DETERMINE  THE  DENSITY  OF  SEAWATER  AT  ATMOSPHERIC  PRESSURE. 

AT  =  AO  +  T*(A1  +  T*(A2  +  T*(A3  +  T*( A4  +  T*(A5  +  T*A6))))) 

BT  =  BO  +  T*(B1  +  T*(B2  +  T*(B3  +  T*B4))) 

CT  =  CO  +  T*  (Cl  +  T*(C2  +  T*C3)) 

DT  =  DO  +  T*(D1  +  T*(D2  +  T*D3)) 

RO  =  AT  +  S* ( BT  +  CT*SQRT(S)  +  DT*S) 

C  DETERMINE  THE  CORRESPONDING  VALUES  OF  SIGMA-T  AND  SPECIFIC  VOLUME. 
SIGMAT  =  1000.0*(RO  -  1.0) 

VO  =  l.O/RO 

RETURN 

END 


V.3  Subroutine  DSW1A 


SUBROUTINE  DSW1A(S,T,K0,A,B) 

C  A  SUBROUTINE  TO  DETERMINE  THE  PRESSURE  CORRECTION  COEFFICIENTS  FOR  THE 
C  SPECIFIC  VOLUME  OF  SEAWATER  AT  SPECIFIED  SALINITY  AND  TEMPERATURE 
C  USING  THE  SOUND-SPEED  DERIVED  EQUATION  OF  STATE  OF  CHEN  AND  MILLERO 
C  (1). 

C  INPUT  PARAMETERS: 

C  S  REAL*4 .  SALINITY  (PARTS/THOUSAND). 

C  T  REAL*4 .  TEMPERATURE  (DEGREES  CENTIGRADE). 

C  OUTPUT  PARAMETERS: 

C  KO ,A,B  REAL*8 .  PRESSURE  CORRECTION  COEFFICIENTS  (SEE  REFERENCE  1). 

C  REFERENCES: 

C  (1)  CHEN,  C.-T.,  AND  MILLERO,  F.J.  1978.  THE  EQUATION  OF  STATE  OF 
C  SEAWATER  DETERMINED  FROM  SOUND  SPEEDS.  J.  MAR.  RES.,  VOL.  36, 

C  657-691. 

C 

REAL*8  KWO.KWl ,KW2 ,KW3 ,KW4 ,KW5 
REAL*8  AWO , AW 1 , AW2 , AW3 , AW4 
REAL*8  BWO , BW1 , BW2 , BW3 , BW4 
REAL*8  AS0.AS1 ,AS2,AS3 
REAL*8  BSO.BSl 
REAL*8  CSO.CSl ,CS2 
REAL*8  DS 

REAL*8  ESO ,ES1 , ES2 

REAL*8  KO , A , B , KW , AW , BW , AS , BS , CS , ES 
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DATA  KWO , KW1 , KW2 ,KWJ , KW4 , KW5/ 1 9652 . 1 7 , 148 . 183 , -2 . 29995 , 1 . 28 ID-02 , 
1  -4 . 9 1564D-05 , 1 . 03553D-07/ 

DATA  AWO , AW1 , AW2 , AW3 , AW4/3 . 26138,5. 223D-04 , 1 . 324D-04 , -7 . 655D-07 , 

1  8.584D-10/ 

DATA  BWO ,  BW  l , BW2 , BW3 , BW4/ 7 . 206 111-05 , -5 . 8948D-06 , 8 . 699D-08 , 

1  -1 .0I0D-09.4.322D-I2/ 

DATA  AS0.AS1 ,AS2 , AS 3/53. 751 , -0 . 4607 , 7 . 030D-03 , -5 . 10711-05/ 

DATA  BSO ,BSl/0 . 2322 , -4 . 838D-03/ 

DATA  CSO ,CS 1 ,CS2/4. 692D-03 , -8 . 387D-05 ,4 . 68D-07/ 

DATA  DS/-1.332D-04/ 

DATA  ESO , ESI ,ES2/- 1 . 41 2D-06 ,9 . 006D-08 , - 1 . 551D-09/ 

C  EVALUATE  THE  PRESSURE  CORRECTION  COEFFICIENTS. 

AS  =  ASO  +  T*(AS1  +  T*(AS2  +  T*AS3)) 

BS  =  BSO  +  T*BS1 

CS  =  CSO  +  T*(CS1  +  T*CS2) 

ES  =  ESO  +  T*(ES1  +  T*ES2) 

KW  =  KWO  +  T*(KW1  +  T*(KW2  +  T*(KW3  +  T*(KW4  +  T*KW5)))) 

AW  =  AWO  +  T*(AW1  +  T*(AW2  +  T*(AW3  +  T*AW4))) 

BW  =  BWO  +  T*(BW1  +  T*(BW2  +  T'-'(BW3  +  T*BW4))) 

S32  =  S*SQRT(S) 

KO  =  KW  +  AS-S  +  BS*S32 

A  =  AW  +  CS* S  +  DS*S32 

B  =  BW  +  ES*S 

RETURN 

END 


V.4  Subroutine CSW2 


SUBROUTINE  DSW2 ( P , VO , KO , A , B , KP , VP , BETA ) 

C  A  SUBROUTINE  TO  DETERMINE  THE  SPECIFIC  VOLUME  AND  ISOTHERMAL 
C  COMPRESSIBILITY  OF  SEAWATER  AT  A  SPECIFIED  PRESSURE,  TEMPERATURE  AND 
C  SALINITY.  IT  IS  ASSUMED  THAT  PREVIOUS  CALLS  TO  DSW1  AND  DSW1A  HAVE 
C  BEEN  MADE, AND  THAT  THE  SPECIFIC  VOLUME  AT  ONE  ATMOSPHERE  PRESSURE, 

C  AND  THE  PRESSURE  CORRECTION  COEFFICIENTS  HAVE  BEEN  COMPUTED. 

C  INPUT  PARAMETERS: 

C  P  REAL*4 .  PRESSURE  (BARS). 

C  VO  REAL*4 .  SPECIFIC  VOLUME  OF  SEAWATER  (CM**3/G)  AT  ONE  ATMOSPHERE 
C  PRESSURE . 

C  KO,A,B  REAL*8 .  PRESSURE  CORRECTION  COEFFICIENTS  (SEE  REFERENCE  1). 

C  OUTPUT  PARAMETERS: 

C  KP  REAL*4 .  SECANT  BULK  MODULUS  (BARS)  AT  SPECIFIED  TEMPERATURE, 

C  PRESSURE  AND  SALINITY. 

C  VP  REAL*4 .  SPECIFIC  VOLUME  OF  SEAWATER  (CM**3/G)  CORRECTED  FOR  THE 
C  EFFECT  OF  PRESSURE. 

C  BETA  REAL*4 .  ISOTHERMAL  COMPRESSIBILITY  OF  SEAWATER  (BAR**-1). 

C  REFERENCE: 

C  (1)  CHEN,  C.-T.,  AND  MILLERO,  F.J.  1978.  THE  EQUATION  OF  STATE  OF 
C  SEAWATER  DETERMINED  FROM  SOUND  SPEEDS.  J.  MAR.  RES.,  VOL.  36, 

C  657-691. 

C 

REAL*4  KP 
REAL*8  KO, A,B 

C  COMPUTE  THE  SECANT  BULK  MODULUS. 
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KP  =  KO  +  P*(A  +  P'-'B) 

C  COMPUTE  THE  SPECIFIC  VOLUME  OF  SEAWATER  AT  ELEVATED  PRESSURE. 
VP  =  VO- (1.0  -  P/KP) 

C  COMPUTE  THE  ISOTHERMAL  COMPRESSIBILITY  OF  SEAWATER. 

BETA  =  VO-  ( KO  -  B*P*P)/(VP*KP*KP) 

RETURN 

END 


V.5  Subroutine  PSW1 


SUBROUTINE  PSW1 (S ,T , P , VO ,VP ,K0 ,A , B .ALPHA ,CP) 

C  A  SUBROUTINE  TO  COMPUTE  THE  EXPANSIBILITY  AND  CONSTANT-PRESSURE 
C  SPECIFIC  HEAT  OF  SEAWATER  AT  SPECIFIED  TEMPERATURE,  PRESSURE  AND 
C  SALINITY.  IT  IS  ASSUMED  THAT  PRIOR  CALLS  HAVE  BEEN  MADE  TO  SUBROUTINES 
C  DSW1 ,  DSW1A  AND  DSW2  TO  COMPUTE  THE  SPECIFIC  VOLUME  OF  SEAWATER  AT 
C  ATMOSPHERIC  PRESSURE  AND  AT  PRESSURE  P,  AND  THE  PRESSURE  CORRECTION 
C  COEFFICIENTS  KO,  A  AND  B. 

C  INPUT  PARAMETERS: 

C  S  REAL-4 .  SALINITY  (PARTS/THOUSAND). 

C  T  REAL-4.  TEMPERATURE  (DEGREES  CENTIGRADE). 

C  P  REAL- 4.  PRESSURE  (BARS). 

VO  REAL "'4 .  SPECIFIC  VOLUME  OF  SEAWATER  (CM--3/G)  AT  SPECIFIED 

SALINITY  AND  TEMPERATURE  AND  ATMOSPHERIC  PRESSURE. 

VP  REAL* 4 .  SPECIFIC  VOLUME  OF  SEAWATER  (CM**3/G)  AT  SPECIFIED 

TEMPERATURE,  PRESSURE  AND  SALINITY. 

KO , A , B  REAL*8 .  PRESSURE  CORRECTION  COEFFICIENTS  (SEE  REKRENCE  1). 
OUTPUT  PARAMETERS: 

ALPHA  REAL*4.  EXPANSIBILITY  OF  SEAWATER  (DEGREES  CENTIGRADE**- 1 )  AT 
SPECIFIED  SALINITY,  TEMPERATURE  AND  PRESSURE. 

CP  REAL*4 .  CONSTANT  PRESSURE  SPECIFIC  HEAT  OF  SEAWATER 

( J . G**- 1 . DEG**- 1 )  AT  SPECIFIED  SALINITY,  TEMPERATURE  AND 
PRESSURE. 

REFERENCE : 

(1)  CHEN,  C.-T.,  AND  MILLERO,  F.J.  1978.  THE  EQUATION  OF  STATE  OF 
SEAWATER  DETERMINED  FROM  SOUND  SPEEDS.  J.  MAR.  RES.,  VOL.  36, 
657-691 . 

REAL*8  DRDT , DVDT , DKDT , DADT , DBDT , X , D2RDT2 , D2VDT2 , D2KDT2 , D2ADT2 , 

1  D2BDT2 ,CPWO , ACP , BCP.CPO ,Q, 110, Ill, 1 20, 121, 1 22, 123, 1 30, 131, 1 32, 

2  133,134, 135, A1,B1,B2,B3, Cl , C2 , C3 , C4 , C5 , I 2 , DSQ , ABP2 , KB2AS , KO , A , B , 

3  TBPA , DSQ , QX , QKO , XSQ , KOSQ , KOD , AD , BD , PD , AERR , RERR , ERROR 
REAL*8  DDER 

C0MM0N/PSW1 A/KOD , AD , BD , DVDT , D2VDT2 , DKDT , D2KDT2 , DADT , D2ADT2 , 

1  DBDT.D2BDT2 
EXTERNAL  DDER 
KOD  =  KO 
AD  =  A 
BD  =  B 
PD  =  P 

AERR  =  1.0D-06 
RERR  =  1.0D-06 
C  EVALUATE  FIRST-ORDER  PARTIAL  DERIVATIVES. 

DRDT  =  6 . 79 14D-05  +  T*(-l .81788D-05  +  T*(3.0513D-07  + 
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t  r'  t -S.  1  18411-09  +  T*(5 .  79GD-  I  I  -  T*3 . 00/5D-  1  ;i ) )  )  )  + 

2  S* (  -4 . 4490D-06  +  T*(2.097D-07  +  T*  (-3 . 774D-09  +  T*1 . 326D- 1 1 )  )  ■+ 

3  SQKT(S) "•'•'( 2 . 844  ID-07  +  T*(-3 . 3742D-08  +  T* 8 . 49774D- 10) )  + 

4  S*(-l  .97975D-08  +  T*(3 . 3282D-09  -  T*9 . 3609D- 1 1 )  )  ) 

DVDT  =  -V0*V0*DRDT 

DKDT  =  1 . 48 183D+02  +  T*(-4.5999D+00  +  T*(3.843D-02  + 

1  T*(-l  .966256D-04  +  T*5 . 1 7765D-07) ) )  +  S*(-4 . 607D-0 1  + 

2  T*( 1 . 406D-02  -  T*1 . 532 ID-04)  -  4 . 838D-03*SQRT(S) ) 

DADT  =  5 . 223D-04  +  T*(2 . 648D-04  +  T*(-2.2965D-06  +  T*3 . 4336D-09) ) 

1  +  S*(-8 . 387D-05  +  T*9.36D-07) 

DBDT  =  -5 . 8948D-06  +  T*( 1 . 7398D-07  +  T*(-3 . 03D-09  + 

1  T*1 .  7288D-  11))  +  S'"  (9 . 006D-08  -  T*3.102D-09) 

C  COMPUTE  EXPANSIBILITY. 

X  =  KO  +  P "  ( A  +  P“B) 

ALPHA  =  DVDT/ VP  -  P'-DVDT/  (VP’-X)  +  P-- VO- (DKDT  +  P* (DADT  + 

1  P"DBDT))/(VP"X"X) 

C  EVALUATE  SECOND-ORDER  PARTIAL  DERIVATIVES. 

D2RDT2  =  -1.81 7880-05  +  T'~(  6.  1026D-07  +  T*(-l -54152D-08  + 

1  T"  (2 . 3 184D- 10  -  T* 1 . 30375D- 12)))  +  S- (2 . 097D-07  +  T"(~7 .548D-09  + 

2  T*3 . 978D- 1 1 )  +  SQRT(S)- (-3 . 3742D-08  +  T* 1 . 699548D-09 )  + 

3  S*(3 . 3282D-09  -  T* 1 .87218D-10) ) 

D2VDT2  =  V j" VO" ( -D2RDT2  +  2 . 0- VO-DRDT-DRDT) 

D2KDT2  =  -4.5999  +  T*( 7 . 686D-02  +  T*(-5 . 898768D-04  + 

1  T--2 . 07  106D-06)  )  +  S* ( 1 . 406D-02  -  T*3 . 0642D-04) 

D2ADT2  =  2 . 648D-04  +  T*(-4. 593D-06  +  T*1 . 03008D-08)  +  S*9.36D-07 
D2BDT2  =  1 .  7398D-07  +  T*(-6 . 06D-09  +  T>^5 . 1864D- 1 1 )  -  S--3.102D-09 
C  COMPUTE  SPECIFIC  HEAT  OF  SEAWATER  AT  ATMOSPHERIC  PRESSURE. 

CPWO  -  4.21 74D+00  +  Pv(-3. 720283D-03  +  Tv(l . 412855D-04  + 

1  T*(-2 - 654387D-06  +  1^-2 . 093236D-08) ) ) 

ACP  =  -7 . 644D-03  +  T* ( 1 .0727D-04  -  P-1.38D-06) 

BCP  =  1.77D-04  +  T‘"  (.-4.08D-06  +  T*5 . 35D-08) 

CPO  =  CPWO  +  S*  (ACP  +  SQRT(S)'-BCP) 

C  EVALUATE  INTEGRAL  IN  THE  SPECIFIC  HEAT  EQUATION.  IF  |B|<0.5D-05, 

C  INTEGRATE  NUMERICALLY,  OTHERWISE  USE  ANALYTICAL  EXPRESSION. 

IF  (DABS(B) .LT.0.5D-05)  GO  TO  1 
Q  =  4 . 0*K0*B  -  A'"A 
TBPA  =  2 . 0-B*P  +  A 
IF  (Q.GT.O.OD+OO)  DSi>  -  DSQRT(Q) 

IF  (Q.LT.O.OD+OO)  USy  =  DSQRT(-Q) 

QX  =  Q'"X 
QKO  =  Q*KO 
XSQ  =  X*X 
KOSQ  =  KO"KO 

IF  (Q.GT.O.OD+OO)  110  =  2 . 0*(DATAN(TBPA/DSQ)  -  DATAN (A/DSQ))/ DSQ 

IF  (Q.LT.O.OD+OO)  110  =  DI,OG( (TBPA  -  DSQ)*(A  +  DSQ)/((TBPA  +  DSQ) 

1  *(A  -  DSQ) ) )/DSQ 

120  =  TBPA/QX  -  A /QKO  +  2.0*B"I10/Q 
Ill  =  0 . 5'"DLOG(X/KO)/B  -  0.5*A*I10/B 

121  =  -0.5»(1.0/X  -  1 . O/KO) /B  -  0 . 5*A*I20/B 

122  =  (-P/X  +  K0*I20)/B 

123  =  (Ill  -  KO* 121  -  A*I22)/B 

130  =  (0 . 5"'(TBPA/XSQ  -  A/KOSQ)  +  3 . 0*B* I20)/Q 

131  =  (-0. 25* ( 1 . O/XSQ  -  l.O/KOSQ)  -  0.5*A*I30)/B 

132  =  (-P/XSQ  -  A* 1 3 1  +  K0*1 30 ) / ( 3 . 0*B) 

133  =  (-0 . 5'"'P'"'P/XSQ  +  KO*  1 3 1 )  /  B 

134  =  (-P**3/XSQ  +  A* 1 33  +  3.0*KO*I32)/B 

135  =  (123  -  KO*I 33  -  A*I34)/B 
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C  EVALUATE  COEFFICIENTS. 

A1  =  -D2VDT2 

B1  =  2 . 0*DVDT*DKDT  +  VO-D2KDT2 
62  =  2 . 0*DVDT*DADT  +  V0*D2ADT2 
B3  =  2 . 0*DVDT*DBDT  +  VO-D2BDT2 
Cl  =  -2 . 0*VO*DKDT*DKDT 
C2  =  -4 . 0 *V 0* D ADT*DKDT 

C3  =  -4 . 0*V0*DBDT*DKDT  -  2 . 0*V0*DADT*DADT 
C4  =  -4.0*V0*DADT*DBDT 
C5  -  -2 . 0*V0*DBDT*DBDT 
C  COMPUTE  SPECIFIC  HEAT. 

12  =  P*D2VDT2  +  A1*I 1 1  +  B1*I21  +  B2*I22  +  B3*I23  +  C1*I31  + 

1  C2*I32  +  C 3*1 33  +  C4*I34  +  C5*I35 
GO  TO  2 

1  CALL  QA05AD(I2, DDER, 0.0D+00, PD, AERR,RERR,1, ERROR, IFLAG) 

IF  (IFLAG.GT.3)  WRITE(6,601)  I FLAG 

601  FORMAT ( ’  POSSIBLE  ERROR  IN  INTEGRAL  EVALUATION,  I FLAG  =  ’,12) 

2  CP  -  CPO  -  0. 1*(T  +  273 . 15)*I2 
RETURN 

END 


V.6  Function  DDER 


DOUBLE  PRECISION  FUNCTION  DDER(P) 

C  A  FUNCTION  TO  EVALUATE  THE  SECOND  DERIVATIVE  OF  SPECIFIC  VOLUME 
C  WITH  RESPECT  TO  TEMPERATURE  AT  A  SPECIFIED  PRESSURE,  P. 

C 

IMPLICIT  REAL*8 ( A-H , 0- Z ) 

COMMON / PSW 1A/K0 , A ,B ,DVDT , D2VDT2 , DKDT , D2KDT2 , DADT , D2ADT2 , 

1  DBDT.D2BDT2 
REAL*8  KO 

X  =  KO  +  P*(A  +  P*B) 

DS1  =  DKDT  +  P1' (DADT  +  P*DBDT ) 

DS2  =  D2KDT2  +  P*(D2ADT2  +  P>D2BDT2) 

DDER  =  D2VDT2  +  (-P-D2VDT2  +  P* ( ( 2 . 0*DVDT*DS 1  +  V0*DS2)  - 
1  2 . 0*V0*DS1*DS1/X)/X)/X 
RETURN 
END 


V.7  Subroutine  SSSW1 


SUBROUTINE  SSSW1 (P,S ,T,U) 

C  A  SUBROUTINE  TO  COMPUTE  THE  SPEED  OF  SOUND  IN  SEAWATER  AT  SPECIFIED 
C  PRESSURE,  SALINITY  AND  TEMPERATURE.  THE  RELATIONSHIP  OF  CHEN  AND 
C  MILLERO  (1)  IS  USED. 

C  INPUT  PARAMETERS: 

C  P  REAL*4 .  PRESSURE  (BARS). 

C  S  REAL-4.  SALINITY  ( PARTS/ THOUSAND ) . 

C  T  REAL*4 .  TEMPERATURE  (DEGREES  CENTIGRADE). 

C  OUTPUT  PARAMETER: 

C  U  REAL*4.  SPEED  OF  SOUND  (M/SEC)  AT  SPECIFIED  PRESSURE,  SALINITY  AND 


C  TEMPERATURE . 

C  REFERENCE: 

C  (1)  CHEN,  C.-T.,  AND  M1LLER0,  F.J.  1977.  SPEED  OF  SOUND  IN  SEAWATER 
C  AT  HIGH  PRESSURES.  J.  ACOUST.  SOC.  AM.,  VOL.  62,  PP.  1129-1135. 

C 

REAL- 8  AHO , AH  1 , AH2 , AH3 , AH4 , AH5 

REAL- 8  BHO , BH 1 , BH2 , BH3 , BH4 

REAL- 8  CHO.CHl ,CH2,CH3,CH4 

REAL- 8  DHO.DH1 ,DH2 

REAL-8  ADO, ADI ,AD2 ,Af)3 ,  AD4 

REA  1.-8  HDO.BDl 

REAL-8  A 1 0  ,  A  1 1  ,A12,.\13,A14 

REAL-8  A20 , A2 1 ,A22 , A23 

REAL* 8  A30 , AJ 1 , A32 

REAL- 8  BO ,B1 

REAL-8  CO 

REAL-8  AH , BH , CH , DH , AD , BD , CD , A 1 , A2 , A3 , A , B , C , UPH20 , UD 

DATA  AHO , AH  1 , AH2 , AH3 , AH4 , AH5/ 1 402. 388,5.03711,-5. 8085 2D- 02 , 

1  3 . 3420D-04 , -  1 . 47800D-06 , 3 . 1464D-09/ 

DATA  BHO.BHl , BH2 , BH3 , BH4/0 . 153563 , 6 . 8982D-04 , -8 . 1788D-06, 

1  1 .36210-07,-6. 1 185D- 10/ 

DATA  CHO , CH 1 , CH2 , CH3 , CH4/3 . 1260D-05 , - 1 . 7 107D-06 , 2 . 5974D-08 , 

1  -2.5335D-10, 1. 0405D- 12/ 

DATA  DHO ,DH1 .DH2/-9 . 7729D-09 , 3 . 8504D- 10 , -2 . 3643D-12/ 

DATA  ADO , ADI , AD2 , AD3 , AD4/ 1.389,-1. 262D-03 , 7 . 164D-05 , 2 . 006D-06 , 

1  -3.21D-08/ 

DATA  BD0,BD1/-1 . 9220-02,-4. 42D-05/ 

DATA  CD/ 1 . 727D-03/ 

DATA  A 10, All ,A12,A13,Al4/9. 4742D-05 , - 1 . 2580D-05 , -6 . 4885D-08 , 

1  1 . 0507D-08,-2.0122D- 10/ 

DATA  A20 ,A21 , A22 , A23/-3 . 9064D-07 ,9 . 1041D-09 , - 1 . 6002D-10 , 7 . 988D- 12/ 
DATA  A30 , A3 1 , A32/ 1 . 100D- 10 , 6 . 649D- 12,-3. 389D- 1 3/ 

DATA  BO , Bl/ 7 . 3637D-05 , 1 . 7945D-07/ 

DATA  C0/-7 . 9836D-06/ 

C  COMPUTE  THE  SPEED  OF  SOUND  IN  PURE  WATER  AT  THE  DESIGNATED  TEMPERATURE 
C  AND  PRESSURE. 

AH  =  AHO  +  T "  ( AH  1  +  T-(AH2  +  T--(AH3  +  T-(AH4  +  T-AH5)))) 

BH  =  BHO  +  T-(BH1  +  T'"'(BH2  +  T*(BH3  +  T-BH4))) 

CH  =  CHO  +  T"'(CH1  +  T"  (CH2  +  T*(CH3  +  T*CH4))) 

DH  =  DHO  +  T"  (DH1  +  T*DH2) 

UPH20  =  AH  +  P*  ( BH  +  P-(CH  +  P*DH)) 

C  COMPUTE  THE  DIFFERENTIAL  BETWEEN  THE  SPEED  OF  SOUND  IN  SEAWATER  AND 
C  IN  PURE  WATER  AT  THE  DESIGNATED  TEMPERATURE  AND  PRESSURE. 

AD  =  ADO  +  T*(AD1  +  T*(AD2  +  T*(AD3  +  T"AD4) ) ) 

BD  =  BDO  +  T*BD 1 

UD  =  S*(AD  +  BD-SQRT(S)  +  CD-S) 

C  COMPUTE  THE  SPEED  OF  SOUND  IN  SEAWATER  AT  DESIGNATED  TEMPERATURE, 

C  PRFSSURE  AND  SALINITY. 

A1  =  AlO  +  r-(All  +  T-(A12  +  T-  ( A 1 3  +  T*Al4))) 

A2  =  A20  +  T*(A2 1  +  T-(A22  +  T*A23)) 

A3  =  A30  +  T*(A31  +  r-A32) 

A  =  P-(A1  +  P"(A2  +  P-A3) ) 

B  =  P-(B0  +  T-Bl) 

C  =  CO-P 

U  =  S ( A  +  B*SQRT(S)  +  C*S)  +  UD  +  UPH20 

RETURN 

END 
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V.8  Function  SAI, 


FUNCTION  SAL(LAT4,LONC4,T4) 

C  A  FUNCTION  TO  COMPUTE  THE  SALINITY  CORRESPONDING  TO  A  GIVEN 
C  TEMPERATURE  IN  A  SPECIFIC  GEOGRAPHICAL  LOCATION  USING  A  STANDARD  T-S 
C  CURVE. 

C  PARAMETERS  PASSED: 

C  I.AT4  REAL-4.  LATITUDE  (DEGREES). 

C  LONG4  REAL- 4 .  LONGITUDE  (DEGREES). 

C  T'  REAL-4 .  TEMPERATURE  (DEGREES  CENTIGRADE). 

C 

IMPLICIT  REAL--8  (A-H.O-Z) 

REAL'"4  LAT 4 ,  LONG4  , T4 

REAL*8  XN(8, 12) ,FN(8, 12) ,GN(8, 12) ,X( 12) ,F( 12) ,G( 12) , LAT, LONG 
INTEGER- 4  NP(8) 

DATA  NP( 1 ) ,NP(2) ,NP(3) ,NP(4) ,NP(5) ,NP(b) ,NP(7) ,NP(8)/ 

I  12,9,12,11,9,5,9,9/ 

DATA  XN ( 1 , 1) ,XN(1 ,2),XN(1 ,3),XN(1 ,4) ,XN( 1 ,5),XN(1 ,6) ,XN(1 ,7) , 

1  XN( 1, 8), XN( 1,9 ), XN( 1, 10), XN ( 1 , 1 1 ) , XN ( 1,12)/ 1.7399998, 3. 4487486, 

2  5. 1574984,6.8662472,8.5749969,10.283745,11.992493,15.409996, 

3  18 . 827484 , 22 . 244995 , 25 . 66249 1,29. 080002/ 

DATA  FN(1,1),FN(1,2),FN(1,3),FN(1,4),FN(1,5),FN(1,6),FN(1,7), 

1  FN(1,8),FN(1,9),FN(1,10),FN(1,11),FN(1,12)/34.726685,34.551025, 

2  34 . 463531 , 34 . 503052 , 34 . 642029 , 34 . 806015 , 34 .975433 , 35 . 352448 , 

3  35 . 620544 , 35 . 566254 , 35 . 252518, 34 . 667648/ 

DATA  GN(1,I),GN(I,2),GN(1,3),GN(1,4),GN(1,5),GN(1,6),GN(1,7), 

1  GN( 1, 8),GN ( 1, 9 ), GN ( 1, 10 ), CN( 1, 1 1 ) , GN ( 1 , 12 )/-0. 11545753, 

2  -0.082896233,-0.014968872,0.058553696,0.094157219,0.096733093, 

3  0.10424614,0.10685921,0.034598351,-0.057549477,-0.12745190, 

4  -0.22149563/ 

DATA  XN(2,1),XN(2,2),XN(2,3),XN(2,4) ,XN(2,5) ,XN(2,6) ,XN(2,7) , 

1  XN (2,8 ) ,XN(2,9)/1 .7799997,5. 1999969,8.6199951 ,12.039993, 

2  15.459999,18.879990,22.299988,25.719986,29.139999/ 

DATA  FN(2,1),FN(2,2),FN(2,3),FN(2,4),FN(2,5),FN(2,6),FN(2,7), 

1  FN( 2,8 ) ,FN(2,9)/34. 699966. 34.441 177,34. 607269,34. 9805 15, 

2  35.330429,35.642120,35.758163,35.500107,34.654861/ 

DATA  GN(2, 1 ) ,GN(2 ,2) ,GN(2 , 3 ) ,GN(2 , 4) , GN(2 ,5 ) , GN(2 ,6) ,GN(2 , 7 ) , 

1  GN(2,8) , GN (2, 9 )/-0. 084634781, -0.02335739 1,0. 096752 167, 

2  0.10944939,0.099792480,0.071731567,-0.011492729,-0.15031433, 

3  -0.35505390/ 

DATA  XN(3,1),XN(3,2),XN(3,3),XN(3,4),XN(3,5),XN(3,6),XN(3,7), 

1  XN(3,8),XN(3,9),XN(3, 10), XN(3, 11) ,XN(3, 12)/ 1.8500004, 3. 4937487, 

2  5.1374979,6.7812462,8.4249954,10.068741,11.712494,14.999992, 

3  18.287476,21.574982,24.862488,28.149994/ 

DATA  FN(3,1),FN(3,2),FN(3,3),FN(3,4),FN(3,5),FN(3,6),FN(3,7), 

1  FN(3,S),FN(3,9),FN(3,10),FN(3,ll),FN(3,12)/34. 682465 ,34.517578, 

2  34 . 406326 , 34 . 443558 ,34.561661,34.729767,34. 925735 , 35 . 267654 , 

3  35.553070,35.674911,35.487518,34.840240/ 

DATA  GN(3,1),GN(3,2),GN(3,3),GN(3,4),GN(3,5),GN(3,6),GN(3,7), 

1  GN(3,8),GN(3,9),GN(3,10),GN(3,ll),GN(3,12)/-0. 090103149 , 

2  -0.097715378,-0.022984505,0.054557800,0.088250160,0.11484337, 

3  0.11683750,0.096542358,0  069462776,-0.0027217865,-0.11839390, 

4  -0.28538132/ 

DATA  XN(4,1),XN(4,2),XN(4,3),XN(4,4),XN(4,5),XN(4,6),XN(4,7), 

1  XN(4,8),XN(4,9),XN(4,10),XN(4,11)/1. 1899996,2.7418737,4.2937479, 
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2  5.84‘><»22  1  ,  / .  (9  /4<)/2  ,1(1.  501244,  I  ).o04<J96, 16. 708740, 19.812500, 

.1  22.91  (>245 , 2h .  020004/ 

DATA  I  N  (4, 1  ) ,  KN(  4,2),  KN (4,1) ,  KN  (4 ,4) ,  FN(4 ,5  ) , KN(4 ,6)  ,FN(4, 7)  , 

1  KN  ( 4 , 8  ) ,  KN  (  4 , 9  ) ,  KN  ( 4  ,  1 0 ) ,  KN  ( 4 , 1  I )/34. 720779,34.626480,34.494614, 

2  34 . 4  72')(>5 , 34 . 59800.', ,  14 . 888580 , 35 . 243  (02 , 35 . 5  1 0086 , 35 . 633240  , 
i  35.389/t,/,  (5.443431/ 

DATA  GN (4 , 1 ) , GN ( 4 , 2 ) , GN ( 4 , 3 ) , GN (4 , 4 ) , GN(4 , 5 ) ,GN (4 ,6) , GN(4 , 7) , 

1  GN  ( 4 , 8  )  ,  GN  ( 4 , 9  .1 ,  GN  (  4 , 1 0  )  ,  GN  ( 4 , 1  1  )  / 0 . 11289215,-0.12562466, 

2  -0.047b 14098, 0.0 18503 189, 0.076808929,0. 12305641 ,0. 10305786, 

3  0.065424919,0.012089/29,-0.036752701,-0.048524857/ 

DATA  XN(5 , 1 ) ,XN (5 , 2 ) , XN(5 , 3 ) ,XN (5 ,4J ,XN(5 , 5 ) ,XN(5 ,6) ,XN(5 , 7) , 

1  XN ( 5 , 8 ) , XN  (.5,91/1. 8999996 ,4.9149990,7. 9299994 , 10 . 944992  , 

2  13.959999, 10.974991 , 19.989990,23.004990,26.020004/ 

DATA  FN ( 5 , 1 ) , KN (5 , 2 1 , KN(5 , 3) ,FN(5 ,4) ,FN(5 ,5) ,FN(5 ,6) ,FN(5 , 7) , 

1  FN(5, 8) ,KHi. 5. 9 1/34.687988 .34. 4729 16 ,34. 588760, 34. 91 1652, 

2  35 . 253723 , 35 . 53082  i ,  (3 . 7004 / 0 , 35 .715652,35. 457870/ 

DATA  GN  (.5,1 )  ,  GN  (.5 , 2  )  ,  GN  ( 5 , 3 )  ,  GN(5 ,4 )  ,  GN  (5 ,5  )  ,  GN(5 ,6)  ,GN(5,7)  , 

1  GN(5,8),GN(5,9)/-0. 072339058, -0.028326035, 0.08690071 1,0. 11 726761, 

2  0. 10508237,0. 076093674,0. 03447 15 12, -0.030061 722,-0. 15561867/ 

DATA  XN ( 6 , 1 ) , XN ( 6 , 2 ) , XN ( 6 , 3 ) ,XN (6 ,4) ,XN(6 , 5 )/ 1 . 8500004 , 7 . 81 74982 , 

1  13.784990, 19.752487,25.720001/ 

DATA  KN ( 6 , 1 ) ,FN(6 ,2) , FN(6 , 3) ,FN(6,4) ,FN(6,5)/34 .691071 ,34.583099, 

1  35.24  7391 ,35.713013,  (5.447403/ 

DATA  GN(6, 1) ,GN(6,2) ,GN (0,3) ,GN(o,4) ,GN(6,5)/-0. 13153839, 

1  0.075304031 ,0. 1 10006  (  S , 0 . 05209  7 182 , -0 . 220237 73/ 

DATA  XN ( 7 , 1) , XN ( 7,2) , XN ( 7 , 3 ) , XN ( 7 , 4 ) , XN ( 7 , 5 ) , XN ( 7 , 6 ) , XN ( 7 , 7 ) , 

1  XN(7,8),XN(7,9)/1. 1049997,3.5224981 ,5.9349966,8.3474884, 

2  10.759995,13. 172493, 15.384999,17.997498,20.410004/ 

DATA  FN ( 7 , 1 ) , FN ( 7 , 2 ) . KN ( 7 ,  5 ) , KN ( 7 , 4 ) , FN ( 7 , 5 ) , FN ( 7 , 6 ) , FN ( 7 , 7 ) , 

1  FN (7, 8 ),FN( 7, 9) /34. 728088,34.538361 ,34.459381 ,34.632278, 

2  34 .937805 , 35 . 23o90S , 35 . 452240 , 35 . 554504 , 35 . 544281/ 

DATA  GN ( 7 , 1 ) , GN( 7 , 2 ) , ON ( 7 , 3) , GN( 7 , 4) , GN ( 7 ,5) , GN( 7 ,6) , GN(7 , 7) , 

1  GN ( 7 , 8 ) ,GN(7,9)/0.044  (59207,-0. 10148716,0.027439117,0.10853386, 

2  0. 13336277,0. 10989761 ,0.066784859,0.017894745,-0.023929596/ 

DATA  XN ( 8 , 1 ) , XN ( 8 , 2 ) , XN ( 8 . 3 ) , XN (8 , 4 ) , XN (8 , 5 ) , XN ( 8 , 6 ) , XN (8 , 7 ) , 

1  XN(8,8) , XN ( 8 , 9 ) / 1 . i 099997 ,4. 0724974 , 7 . 034995 1 ,9 . 9974899 , 

2  12 . 95999 1,15. 9224°  '  18.884995,21.847488,24.809998/ 

DATA  FN(8 , 1 ) , FN (8 , .  N  .  , 3) , FN (8 , 4) , FN (8 ,5 ) ,FN(8 ,6) ,FN(8 , 7) , 

1  FN(8 , 8) , FN (8 , 9 ) /34 . /420  (5 , 34 . 507553 , 34 . 523254 ,34 . 843521 , 

2  35 . 205383 , 35 . 45  3888  ,  33 .  580 109 ,35.633163,35. 646423/ 

DATA  GN ( 8 , 1 ) , GN ( 8 , 2 ) , GN ( 8 , 3 ) , GN ( 8 , 4 ) ,GN(8,5) ,GN(8,6) ,GN(8,7) , 

1  GN ( 8 , 8 ) ,GN(8, 9 )/-(). 02 7276993, -0.0659 16061 ,0.069396019, 

2  0. 12857723,0. 10705471 ,0.061277390,0.027302742,0.011056900, 

3  -0.0043/54578/ 

LAT  =  DB1.KUAT4) 

LONG  -  DBLE(L0NG4) 

T  =  DBLE(T4) 

IF  (LAT. LE. 25. .AND. LONG. LE . 160.)  J  =  1 

IF  (LAT. I.E. 25. .AND. LONG. GT. 160. .AND. LONG. LE. 170.)  J  =  2 

IF  (LAT. LE. 25. .AND. LONG. GT. 170. )  J  -  3 

IF  (LAT. GT. 25. .AND. LAT. LE. 35. .AND. LONG. LE. 160.)  J=4 

IF  (LAT.GT.25. . AND . LAT . LE . 35 . . AND . LONG . GT. 160 . . AND. LONG . LE . 1 70 . ) 

1  J  =  5 

IF  (LAT.GT.25. .AND. LAT. LE. 35. .AND. LONG.GT. 170.)  J  =  6 
IF  (LAT. GT. 35. .AND. LONG. LE. 160.)  J  =  7 
IF  (LAT.GT. 35. .AND. LONG.GT. 160. .AND. LONG. LE. 170. )  J  =  8 
N  =  NP(J) 
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DO  1  1  =  1,  N 
X(l)  =  XN(J,I) 

F(I)  =  FN(J.I) 

1  G(  I )  =  GN(J, I) 

IF  (T.GE.X(l).AND.T.LE.X(N))  GO  TO  2 
IF  (T.LT.X(l))  SAL  =  F(l)  +  (T  -  X(1))*G(1) 
IF  (T.GT.X(N))  SAL  =  F(N)  +  (T  -  X(N))*G(N) 
RETURN 

2  SAL  =  TG01BD(-1 ,N,X,F,G,T) 

RETURN 

END 


V.9  Subroutine  PLOT1 


SUBROUTINE  PLOT 1 (DEPTH , TEMP , SAL , S IGT , VELOC , L , I J , ISCALE ) 

DIMENSION  DEPTH (200) ,  TEMP(200) ,SAL(200) ,SIGT(200) ,VEL0C(200) , 

1  T(200) ,D(200) ,S(200) ,V(200) ,SG(200) ,Z1 (5) ,Z2(5) ,ZP(5) 

DATA  Zl/’  500','  375','  250','  125','  0  '/, 

1  Z2/ ' 1000 ' , '  750’,'  500’,’  250','  0  '/ 

CALL  PLOTS (IBUF, NO, NUM) 

IF( ISCALE. EQ. 2)  GO  TO  20 

DO  10  1=1,5 
10  ZP(I)  =  Z1(I) 

GO  TO  40 


20  DO  30  1=1,5 

30  ZP(I)  =  Z2(I) 

40  CONTINUE 

DEPTH(L+1 )=0. 

DEPTH(L+2)=  125.*  FLOAT (I SCALE) 
TEMP(L+1)=0. 

TEMP(L+2)=10. 

SAL(L+1)=34.0 
SAL(L+2)=0. 6666667 
SIGT(L+1 )=23 . 

SIGT(L+2)=1. 3333333 
VELOC(L+l )=1480. 

VE  LOC ( L+2 ) =26 . 666666  7 
CALL  PL0T(0 . ,8.25,22) 

CALL  PLOT(ll .65,0. ,22) 

CALL  PLOT(0 . ,-8.25,22) 

CALL  PLOT (-11.65,0. ,22) 

CALL  PLOT  (1., 1.0, 23) 

K=1 

CALL  PL0T(0. ,3. ,2) 

CALL  PLOT(-0.05,3. ,2) 

CALL  SYMB0L(-0. 1,2.9, . 105,2H30,90. ,2) 
CALL  PL0T(0. ,2. ,3) 

CALL  PLOT (-.05,2. ,2) 

CALL  SYMBOL (-0. 1,1.9,. 105 ,2H20,90. ,2) 
CALL  PLOT(0. , 1 . ,3) 

CALL  PLOT (-0. 05 , 1 . ,2) 

CALL  SYMBOL(-0 . 1 ,0.9, . 105,2H10,90. ,2) 
CALL  PLOT  (0. ,0. ,3) 


CAM.  SYMBOL(-0.3,0. 72,  .  105 , 15HTEMPKKA TURK  (C),90.  ,15) 

CALL  PL0T(0. ,0. ,3) 

50  CALL  PLOT(4 . ,0. ,2) 

CALL  PLOT(4. ,-0.05,2) 

CALL  SYMBOL(4. 20,-0. 1 , . 105 ,ZP( 1 ) , 180 . ,4) 

CALL  PLOT (3 . ,0. ,3) 

CALL  PLOT ( 3 . ,-0.05,2) 

CALL  SYMBOL (3. 20,-0. 1, . 105 , ZP(2) , 180 . ,4) 

CALL  PLOT(2 . ,0. ,3) 

CALL  PLOT (2 . ,-0.05,2) 

CALL  SYMBOL (2 .15,-0.1,.105,ZP(3),180. ,4) 

CALL  PLOT ( 1 . ,0. ,3) 

CALL  PLOT ( 1 . ,-0.05,2) 

CALL  SYMBOL( 1.15, -0.1,. 105, ZP(4), 180. ,4) 

CALL  PLOT(0. ,0. ,3) 

CALL  PLOT (0 . ,-0.05,2) 

CALL  SYMBOL (+0 .05,-0. 1 , . 105 ,ZP(5) , 180 . ,4) 

CALL  SYMBOL(2. 47,-0. 30, . 105, 9HDEPTH  (M),180.,9) 

GO  TO  (60, 70, 80, 90), K 

60  CALL  LINE(DEPTH, TEMP, L, 1,0,0) 

CALL  PLOT  (0.  ,3.9,23) 

CALL  PLOT  (0. ,3. ,2) 

CALL  PLOT  (-0.05,3. ,2) 

CALL  SYMB0L(-0. 1,2.9, . 105,2H36,90. ,2) 

CALL  PL0T(0. ,2.25,3) 

CALL  PL0T(-0 .05,2.25,2) 

CALL  SYMBOL (-0 .1,2.15,.! 05 ,4H35 .5 ,90 . ,4) 

CALL  PLOT(0. ,1.5,3) 

CALL  PLOT(-0. 05, 1.5,2) 


CALL  SYMBOL (-0 . 1 , 1 . 4 , . 105 , 2H35 ,90 . , 2 ) 

CALL  PL0T(0 . ,0.75,3) 

CALL  PL0T(-0. 05, 0.75,2) 

CALL  SYMBOL(-0. 1,0.65, . 105 ,4H34 . 5 ,90 . ,4) 

CALL  PLOT(0. ,0. ,3) 

CALL  PLOT  (-0.05,0. ,2) 

CALL  SYMBOL(-0. 1,-0. 10, . 105,2H34,90. ,2) 

CALL  SYMBOL(-0 .3, 0.77, .105, 14HSALINITY  (PPT) ,90. , 14) 

CALL  PLOT  (0. ,0. ,3) 

K=K+1 

GO  TO  50 

70  CONTINUE 

CALL  LINE  (DEPTH, SAL, L, 1,0,0) 

CALL  PLOT(5. ,0. ,23) 

CALL  PL0T(0. ,3. ,2) 

CALL  PL0T(-0.05,3. ,2) 

CALL  SYMBOL (-0 .1,2.9,. 105 ,2H27 ,90. ,2) 

CALL  PLOT(0. ,2.25,3) 

CALL  PLOT (-0.05,2.25,2) 

CALL  SYMBOL(-0. 1,2.15,. 105,2H26,90. ,2) 

CALL  PL0T(0. ,1.5,3) 

CALL  PLOT (-0.05,1.5,2) 

CALL  SYMB0L(-0. 1,1.4,. 105, 2H25, 90. ,2) 

CALL  PL0T(0. ,0.75,3) 

CALL  PL0T(-0. 05, 0.75,2) 

CALL  SYMBOL(-0. 1 ,0.65 , . 105 ,2H24,90. ,2) 
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CALL  PL0T(0. ,0. ,3) 

CALL  PLOT(- .05,0. ,2) 

CALL  SYMB0L(-0. 1 ,-0.05, . 105,2H23,90. ,2) 

CALL  SYMB01,(-0 .3,1 . 13, . 105 , 7HSIGMA-T.90 . , 7) 

CALL  PLOTtO. ,0. ,23) 

K-K+l 

GO  TO  50 

80  CONTINUE 

CALL  LINE (DEPTH, SIGT,L, 1,0,0) 

CALL  PLOT(0. ,-3.9,23) 

CALL  PLOT (0 . ,3. ,2) 

CALL  PLOT(-0.05 ,3. ,2) 

CALL  SYMBOL(-0. 1,2.8, . 105 ,4H1560 ,90 . ,4) 

CALL  PLOT(0. ,2.25,3) 

CALL  PLOT(-0 .05,2.25,2) 

CALL  SYMBOL (-0 .1,2.15,. 105 ,4H1540 ,90 . ,4) 

CALL  PL0T(0 . ,1.5,3) 

CALL  PLOT(-0. 05, 1.5,2) 

CALL  SYMB0L(-0. 1,1.3,. 105, 4H1520, 90. ,4) 

CALL  PLOT(0. ,0.75,3) 

CALL  PLOT (-0.05 ,0.75,2) 

CALL  SYMB0L(-0. 1,0.55, .105, 4H1500, 90. ,4) 

CALL  PL0T(0 . ,0. ,3) 

CALL  PLOT( -0.05,0. ,2) 

CALL  SYMBOL (-0. 1 ,-0.20, . 105 ,4H1480,90. ,4) 

CALL  SYMBOL (-0.3,0. 5, . 105, 20HS0UND  VELOCITY  (M/S) ,90. ,20) 

CALL  PL0T(0. ,0. ,3) 

K  =  K+l 

GO  TO  50 

90  CALL  LINE (DEPTH, VELOC.L, 1,0,0) 

CALL  SYMBOL (5 . ,0 . 75 , . 21 , 17HC0NSECUTIVE  NO.  ,90. ,17) 

A I J=DFLOAT ( I J ) 

CALL  NUMBER (5 . ,999 . , . 2 1 ,AIJ ,90 . ,-l) 

CALL  ENDPLT 

RETURN 

END 


V.10  Subroutine  LNGPLT  (additional  ENTRY  ADDPLT ) 


S  UBROUT I NE  LNGPLT ( S V , UP I , XT I TLE , XARRAY , YARRA Y , L , I J , I CALL , MANY ) 
REAL*4  XTITLE(5) ,YTITLE(4) ,XARRAY( 150) , YARRAY(150) 

DATA  YTITLE/'DEPT' ,'H  -  ','METR’,'ES  '/ 

IF  (ICALL.EQ. 1)  CALL  PLOTS(IBUF,NO,NUM) 

C 

C  DESIGNATES  ORIGIN 
C 

CALL  PLOT(  1.  ,9.  ,23) 

C 

C 

C  PLOT  VERTICAL  SCALE 
C 

CALL  PLOT(-0. 1  ,  0.  ,2) 

CALL  PLOT(  0.  ,0.  ,2) 
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CALI.  I*L0T(  0.  ,-0.5  ,2) 

CALL  PI.0T(  0.  ,-0.5  ,2) 

X  =  0. 

Y  =  -1. 

DO  10  I  =  1,8 
CALL  PL0T(X,Y,2) 

X  =  -0.1 

CALL  PL0T(X,Y,2) 

X  =  0 

CALL  PL0T(X,Y,2) 

Y  =  Y-0.5 

CALL  PL0T(X,Y,2) 

X  =  -0.05 
CALL  PL0T(X,Y,2) 

X  =  0 

CALL  PL0T(X,Y,2) 

Y  =  Y-0.5 
IQ  CONTINUE 

CALL  PL0T(X,Y,2) 

X  =  -0.1 

CALL  PL0T(X,Y,2) 

PRINT  VERTICAL  TITLE 

CALL  SYMBOL(-0.85 ,-5 .48,0. 14 , YTITLE ,90 . , 16) 

PLOT  NUMBERS  ON  VERTICAL  SCALE. 

PLOT  NUMBERS  SUB-ROUTINE  CALLS  'PLOT  SYMBOL'  AND  REFERS  TO 
X  &  Y  ERRORS  AS  BEING  IN  PLOT  SYMBOL. 

FPN  =  0. 

CALL  NUMBER (-0. 34, -0 .07,0  14,F».  N,0. ,-l) 

I FPN  =  FPN+100 . 

RE-DEFINE  Y  AS  ZERO 

Y  =  0 

DO  20  I  =1,9 

Y  =  Y-l 
FPN  =  I FPN 

CALL  NUMBER (-0. 62, Y,0. 14,FPN,0. ,-l) 

IFPN  =  IFPN+100 
20  CONTINUE 

PLOT  HORIZONTAL  SCALE 

CALL  PLOT(0. ,0.  ,3) 

CALL  PL0T(0. ,0.1,2) 

CALL  PLOT ( 0 . , 0 .  ,2) 

X  =  0. 

Y  =  0. 

DO  30  I  =  1,3 
X  =  X+1.0 
CALL  PLOT(X,Y,2) 

Y  =  Y+0.1 

CALL  PLOT(X,Y,2) 

Y  =  Y-0.1 
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JO  CALL  P1.0T(X,Y,2) 

PRINT  HORIZONTAL  TITLE 
CALL  PL0T(0. ,0. ,3) 

CALL  SYMBOL(0 . 1 ,0. 75 ,0. 14,XTITLE,0. ,20) 

PLOT  NUMBERS  ON  HORIZONTAL  SCALE.  ONLY  THE  START  VALUE 
AND  THE  NEXT  THREE  VALUES  ARE  REQUIRED. 

ISV  =  SV 
IUP1  =  UP1 
IFPN  =  ISV-IUPI 
DO  40  I  =  1,4 
XI  =  1-1 

IFPN  =  IFPN+IUPI 

IF(IFPN.GE.O.AND. IFPN.LT. 10)  XM  =  XI-0.07 
IFCIFPN.GE. 10.AND.IFPN.LT. 100)  XM  =  XI-0.14 
IF(IFPN.GE. 100.AND.IFPN.LT. 1000)  XM  =  XI-0.21 
IF(IFPN.GE. 1000.AND.IFPN.LT. 10000)  XM  =  XI-0.28 
FPN  =  IFPN 

40  CALL  NUMBER(XM, 0.2, 0.14, FPN, 0.,-l) 

MOVE  X  ORIGIN  PROGRESSIVELY  TO  THE  RIGHT  AND  EXTEND  X  GRID 
ENTRY  ADDPLT(XARRAY,YARRAY,L, IJ, ICALL) 

IF(ICALL.EQ.l)  GO  TO  50 

CALL  PLOT(0.5,0.0,23) 

CALL  PLOT ( 2 . 5 , 0 . 0 ,  3) 

CALL  PL0T(3. 0,0.0,  2) 

CALL  PLOT ( 3 . 0 , 0 . 1 ,  2) 

TO  PLOT  EACH  TRACE. 

50  XARRAYCL+l)  =  SV 
XARRAY(L+2)  =  UPI 
YARRAYCL+l)  =  0. 

YARRAYCL+2)  =  -100. 

CALL  LINE(XARRAY,YARRAY,L, 1,0,0) 

TO  PLOT  A  SERIAL  NO.  AT  THE  END  OF  EACH  TRACE 

FIJ  =  IJ 

X  =  (XARRAY(L)-SV)/UPI 
Y  =  YARRAY(L)/(-100. )-0.24 
IFCFIJ.GE.O. AND.FIJ.LT. 10)  X  =  X-0.07 
IFCFIJ.GE.IO.AND.FIJ.LT. 100)  X  =  X-0.14 
IFCFIJ.GE.100.AND.FIJ.lt. 1000)  X  =  X-0.21 
CALL  NUMBERCX.Y.O. 14, FIJ, 0. ,-l) 

IF  C ICALL. EQ. MANY)  CALL  ENDPLT 

RETURN 

END 
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APPENDIX  VII  SAMPLE  CARD  IMAGE  INPUT  FILE 


03 

WSRL  CRUISE  TASMAN  SEA  MP1/79/D1 

08/03/79  1400  151.92  DEG  E  33.72  DEG  S 

657. M  WET  17.8  DRY  22.8  DIR  190  12  KNOTS  1018.1MB 
1  76 

0.0  25.20  5.00  25.18  10.00  25.16  15.00  24.75  20.00  24.13  25.00  23.43 

30.00  22.67  35.00  21.86  40.00  21.14  45.00  20.33  50.00  19.77  55.00  19.38 

60.00  19.00  65.00  18.44  70.00  17.88  75.00  17.39  80.00  17.33  85.00  17.10 

90.00  16.71  95.00  16. 36100. 00  16.00110.00  15.44120.00  14.93130.00  14.57 

140.00  14.21150.00  13.93160.00  13.74170.00  13.56180.00  13.37190.00  13.19 

200.00  13.00210.00  12.86220.00  12.72230.00  12.58240.00  12.44250.00  12.30 

260.00  12.17270.00  12.04280.00  11.86290.00  11.65300.00  11.45310.00  11.24 

320.00  11.04330.00  10.91340.00  10.81350.00  10.70360.00  10.38370.00  10.06 


380.00 

9.94390.00 

9 . 86400 . 00 

9.78410.00 

9.71420.00 

9.63430.00 

9.55 

440.00 

9.48450.00 

9.40460.00 

9.27470.00 

9.14480.00 

9.01490.00 

8.93 

500.00 

8.85510.00 

8.78520.00 

8.70530.00 

8.63540.00 

8.55550.00 

8.47 

560.00 

8.40570.00 

8.32580.00 

8.24590.00 

8.17600.00 

8.09610.00 

8.02 

620.00 

7.69630.00 

7.31640.00 

6.96650.00 

6.76 

WSRL  CRUISE  TASMAN  SEA  MP1/79/D1 
8/03/79  1500  152.17  DEG  E  33.68  DEG  S 

1808. M  WET  17.8  DRY  22.8  DIR  170  14  KNOTS  1018.1MB 

2  99 

0.0  25.80  5.00  25.79  10.00  25.78  15.00  25.77  20.00  25.76  25.00  25.75 

30.00  25.74  35.00  25.73  40.00  25.72  45.00  25.71  50.00  25.70  55.00  24.33 

60.00  22.78  65.00  22.22  70.00  21.70  75.00  21.62  80.00  21.53  85.00  21.33 

90.00  21.06  95.00  20.76100.00  20.47110.00  19.83120.00  19.00130.00  17.79 

140.00  17.07150.00  16. 18160.00  15.68170.00  15.28180.00  14.75190. 00  13.97 

200.00  13.72210.00  13.47220.00  13.22230.00  12.99240.00  12.90250.00  12.80 

260.00  12.42270.00  12.04280.00  11.83290.00  11.63300.00  11.44310.00  11.25 

320.00  11.06330.00  10.90340.00  10.75350.00  10.60360.00  10.41370.00  10.22 

380.00  10.04390.00  9.84400.00  9.64410.00  9.44420.00  9.24430.00  9.04 

440.00  8.87450.00  8.70460.00  8.61470.00  8.52480.00  8.42490.00  8.33 

500.00  8.24510.00  8.15520.00  8.06530.00  7.97540.00  7.88550.00  7.80 

560.00  7.72570.00  7.64580.00  7.55590.00  7.47600.00  7.39610.00  7.31 

620.00  7.22630.00  7.14640.00  7.06650.00  6.98660.00  6.90670.00  6.81 

680.00  6.73690.00  6.65700.00  6.57710.00  6.49720.00  6.41730.00  6.33 

740.00  6.25750.00  6.17760.00  6.09770.00  6.01780.00  5.94790.00  5.88 

800.00  5.81810.00  5.75820.00  5.69830.00  5.62840.00  5.56850.00  5.49 

860.00  5.43870.00  5.36880.00  5.30 

WSRL  CRUISE  TASMAN  SEA  MP1/79/D1 
8/03/79  1600  152.36  DEG  E  33.64  DEG  S 

2848. M  WET  17.8  DRY  22.8  DIR  150  12  KNOTS  1016.3MB 

3  96 

0.0  25.60  5.00  25.55  10.00  25.51  15.00  25.46  20.00  25.41  25.00  25.37 

30.00  25.32  35.00  24.33  40.00  23.75  45.00  23.44  50.00  23.13  55.00  22.86 

60.00  22.64  65.00  22.41  70.00  22.18  75.00  21.96  80.00  21.76  85.00  21.56 

90.00  21.36  95.00  21.16100.00  20. 96110. 00  20.54120.00  20.13130.00  19.63 

140.00  19.11150.00  18.56160.00  18.00170.00  17.55180.00  17.09190.00  16.78 

200.00  16.74210.00  16.59220.00  16.05230.00  15.67240.00  15.30250.00  14.90 

260.00  14.63270.00  14.35280.00  14.08290.00  13.75300.00  13.39310.00  13.04 

320.00  12.79330.00  12.56340.00  12.33350.00  12.10360.00  11.94370.00  11.76 

380.00  11.57390.00  11.39400.00  11.20410.00  11.02420.00  10.71430.00  10.39 

440.00  10.06450.00  10.00460.00  9.81470.00  9.62480.00  9.43490.00  9.25 
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500.00 

9.06510.00 

8.93520.00 

8.82530.00 

8.71540.00 

8.61550.00 

8.50 

560.00 

8.39570.00 

8.29580.00 

8.18590.00 

8.07600.00 

7.97610.00 

7.88 

620.00 

7.78630.00 

7.69640.00 

7.59650.00 

7.50660.00 

7.41670.00 

7.31 

680.00 

7.22690.00 

7.12700.00 

7.03710.00 

6.96720.00 

6.89730.00 

6.83 

740.00 

6.76750.00 

6.70760.00 

6.63770.00 

6.57780.00 

6.51790.00 

6.44 

800.00 

6.38810.00 

6.31820.00 

6.25830.00 

6.19840.00 

6.12850.00 

6.06 

f 

i 

i 


Li 
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APPENDIX  VIII  INPUT  DATA  FORMAT  FOR  AN  XBT  DATA  POINT  RECORD  AS 

USED  BY  AODC 

VIII.  1  Input  record  format 


Element  Name 
and  Level 

Position 
Byte  Nos 

Number 

Bytes 

Usage  and  Meaning 
of  Element 

File  I.D. 

t 

1-2 

2 

Char(2) 

File  I.D.  Always  17 

QUADRT1 

3 

1 

Char(l) 

WMO  3333 

TEN  SQ1 

3 

Char (3) 

Ten  Degree  Square  (WMO) 

FIVE  SQ1 

7 

1 

Char(l) 

Five  Degree  Square 

TWO  SQ1 

8-9 

2 

Char (2) 

Two  Degree  Square 

ONE  SQ1 

10-11 

2 

Char (2) 

One  Degree  Square 

DATE 

YEAR 

12-15 

4 

Char(4) 

Prefix  '19' 

MONTH 

16-17 

2 

Char (2) 

01-12 

DAY 

18-19 

2 

Char (2) 

01-31 

TIME 

HOUR 

20-21 

2 

Char(2) 

00-23 

MIN 

22-23 

2 

Char (2) 

00-59 

REF.  I.D. 

CNTRY 

24-26 

3 

Char(3) 

NODC  CODES 

REF.  NO 

27-31 

5 

Char (5) 

CONSEC 

32-34 

3 

Char(3) 

SHIP 

35-40 

6 

Char(6) 

LATITUDE 

LAT.  DEG 

41-42 

2 

Char(2) 

00-90 

LAT.  MIN 

43-44 

2 

Char (2) 

00-59 

LAT.  HEM 

45 

1 

Char(l) 

N  or  S 

LONGITUDE 

LONG  DEG 

46-48 

3 

Char(3) 

000-179 
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LONG  MIN 

49-50 

2 

Char(2) 

LONG  HEM 

51 

1 

Char(l ) 

NODC  CODES 

BOTTOM 

52 

1 

CharO  ) 

DIGMTH2 

53-54 

2 

Char(2) 

INTER2 

55-56 

2 

Char(2) 

TRESTO2 

57-58 

2 

Char(2) 

OPERATOR 

1NIT 

59-61 

3 

Char (3) 

TRIAL 

62 

1 

Char(l) 

CALDEP 

63-65 

3 

Char (3) 

CALTEM 

66-68 

3 

Char (3) 

INSTRUMENT 

69 

1 

Char(l) 

GRID/33 

70 

1 

Chard ) 

ORIG  CR  NO 

71-78 

8 

Char (8) 

DNP 

79 

1 

Char ( 1 ) 

SKIP 

80 

1 

Char ( 1 ) 

LENGTH 

81-84 

4 

Char (4) 

DEPTH  &  TEMP 

SURTEM 

85-88 

4 

Char(4) 

DEPTH  1 

89-92 

4 

Char(4) 

TEMP  1 

93-96 

4 

Char(4) 

DEPTH  2 

97-100 

4 

Char (4) 

TEMP  2 

101 

4 

Char (4) 

00-59 
E  or  W 

B  =  probe  hit  bottom,  blank  =  did  not 

Method  of  digitisation 

Interval  of  digitisation 

Method  of  treatment  &  stors.,:  of 
initial  points. 

Operator's  Initials 

Number  of  attempts  at  digitising 
trace . 

Depth  of  calibration  tick  in 
Units  of  analogue  grid. 

Negative  depth  value  for 

calibration  tick  above  surface, 
Positive  depth  value  for 

calibration  tick  below  surface. 

Temperature  of  calibration  tick 
in  units  of  analogue  grid. 

16.7°  is  nil-calibration  value. 

1=XBT ,  2=HXBT,  3=SXBT,  4=AXBT 

Grid  codes  of  instrument  1-9. 

Originators  cruise  number. 

Declared  National  Program 

Blank 

Number  of  temperature  values. 

Temperature  at  zero  depth. 

First  depth  to  whole  metres. 

Temp  at  first  depth. 

Second  depth. 

Temperature  at  second  depth. 
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DM  I’ I'll  (N)  4  Ch.ir(4)  Last  depth  Position  =  ( N  —  1  )  +  8r> 

TEMI’  (N)  4  Char (4)  Last  temperature 

Position  =  ( N - 1 )  +  89 


Notes : 

1.  The  location  of  an  XBT  station  on  the  Earth's  surface  is  identified  to 
within  a  10°  square  using  the  WMO  Code  3333  quadrant  and 

10°  square  numbering  system.  Under  this  system,  position  is 
designated  by  a  four-digit  number;  the  first  digit  identifies  the 
quadrant  of  the  globe  (figure  8),  and  may  take  values  of  1(NE),  3(SE), 
5(SW)  or  7(NW);  the  second  digit  is  the  tens  digit  of  the  latitude  of 
the  station;  the  third  and  fourth  digits  are  the  hundreds  and  tens 
digits  of  degrees  of  longitude.  Position  within  a  10°  square  is 
identified  using  the  modified  Canadian  subdivision  systems;  subdivisions 
of  5-,  2-,  and  1-degree  squares  are  coded  as  shown  in  figure  9. 

2.  See  Section  VII  1. 3  for  codes. 

3.  See  Section  VIII. 2  for  codes. 


VI 1 1. 2  BT  instrument  and  grid  code 


Instrument 
and  grid 

Instrument 

code 

Grid 

code 

Shipboard  XBT 

1 

2500  foot  (Ft  &  F) 

1 

750  metre  (M  &  C) 

2 

1500  foot  (Ft  &  F) 

3 

450  metre  (M  &  C) 

4 

1830  metre  (M  &  C) 

5 

Helicopter  XBT 

2 

Unknown 

6 

Submarine  XBT 

3 

Unknown 

7 

Airborne  XBT  (AXBT) 

4 

1000  (Ft  &  F) 

8 
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VIII. 3  Digitisation  method,  interval  of  digitisation,  and 
Data  treatment  and  storage  codes 


I  DIGITISATION  METHOD  Code 

a.  Manual  01 

b.  A-D  conversion  from  original  02 

c.  A-D  conversion  from  copies  03 

d.  Optical  scanning  04 

e.  Direct  digital  output  05 

II  INTERVAL  OF  DIGITISATION  Code 

A.  Fixed  interval 

a.  Sim;  and  S0.1°C  01 

b.  >lm  but  S3m;  and  S0.1°C  02 

c.  >3m  but  S6m;  and  S0.1°C  03 

d.  >6m;  and  S0.1°C  04 

e.  Sim;  and  S0.2°C  11 

f.  >lm  but  S3m;  and  S0.2°C  12 

g.  >3m  but  S6m;  and  SO. 2 °C  13 

B.  Variable  (flexure  points) 

a.  Manually  determined  31 

b.  Statistically  determined  32 

c.  Physically  determined  33 

C.  Combination  of  fixed  and  variable 

a.  Every  3  ft  to  900  ft  then  flexure  points  below  900  ft 

III  TREATMENT  AND  STORAGE  OF  DIGITISED  DATA  Code 

A.  Single  digitisation. 

a.  No  treatment,  stored  as  digitised  01 

DATA  COMPRESSION  resulting  in 

b.  Fit  within  0.05°C  02 

c.  Fit  within  0.1°C  03 

d.  Fit  within  0.2°C  04 

e.  Fit  within  0.3°C  05 

f.  Fit  within  0 . 7°C  06 

B.  Dual  digitisation  and  averaging 

a.  No  treatment,  stored  as  digitised  21 

DATA  COMPRESSION  after  averaging 

b.  Fit  within  0 . 05°C  22 

c.  Fit  within  0.1°C  23 

d.  Fit  within  0.2°C  24 

e.  Fit  within  0.3°C  25 

f.  Fit  within  0.5°C  26 


C.  Data  points  at  fixed  intervals  or  selected  intervals  retained 

and  stored 


27 
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APPENDIX  IX  HARWELL  SUBROUTINES  USED  BY  THE  DATA  REDUCTION  PROGRAMME 

The  XBT  data  reduction  programme  uses  various  subroutines  from  the  Harwell 
Subroutine  Library  for  spline  fitting,  and  for  interpolation  and  quadrature 
using  cubic  spline  fits.  The  following  excerpts  from  the  Harwell  Subroutine 
Library  User's  Manual  describe  the  functions  and  argument  lists  for  the 
subroutines  used. 


IX. 1  Subroutine  TB04AD 


1 .  Purpose 

Given  function  values  f1)f2,...,f  at  points  X!<x2<. . .<x  which  need  not  be 
equally  spaced,  this  routine  finRs  a  cubic  spline  S(x)n  that  interpolates 
the  function  values,  ie 


S(x.) 


,  i=l  ,2, .  .  .  ,n 


where  S(x) 
of- freedom 
derivat ive 
defined  on 
knots,  and 


has  knots  at  the  points  Xx,X2,...,x  .  The  extra  two  degrees- 
for  such  an  interpolation  are  taken  up  by  forcing  the  third 
of  S(x)  at  the  points  x  and  x  ^  to  be  continuous.  S(x)  is 
return  from  the  routine  By  n,  tRe  knots  x.,  its  values  f.  at  the 
the  values  d.  of  its  first  derivatives  at1the  knots. 

t 


2.  Argument  list 


SUBROUTINE  TB04AD(N,X,F,D, W) 


The  arguments  N,  X  and  F  must  be  set  by  the  user  before  calling  TB04AD. 

N  is  INTEGER*4  and  should  be  set  to  n,  the  number  of  points. 

X  is  a  REAL*8  array  of  length  at  least  N  containing  the  points  x., 

i=l,...,n  (which  are  also  the  knots).  These  must  be  ordered  such  that 
x!<x2<...<x  .  If  this  condition  is  not  fulfilled,  W(l) 
is  set  to  1  and  a  diagnostic  message  is  printed  by  the  routine. 

F  is  a  REAL*8  array  of  length  at  least  N  which  must  contain  the  function 
values,  fj,f2,...,fn. 

D  is  a  REAL*8  array  of  length  at  least  N  in  which  the  routine  puts  the 
values  of  the  first  derivatives  of  the  spline  S(x)  at  the  knots  x^. 

W  is  a  REAL*8  array  of  length  at  least  3*N  used  by  the  routine  as  work¬ 
space.  W(l)  is  set  to  zero  on  a  successful  return  and  to  1  on  a 
failure. 


* 
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IX. 2  Function  TG01BD 


1 .  Purpose 

Given  a  cubic  spline  in  terms  of  its  knots,  values  and  first  derivative  at 
the  knots,  this  function  calculates  the  value  of  the  spline  at  any  other 
point . 

If  the  point  at  which  the  value  is  required  lies  outside  the  interval 
between  the  first  and  last  knots,  the  value  of  the  spline  is  set  to  zero. 

2.  Argument  list 

DOUBLE  PRECISION  FUNCTION  TG01BD(I ,N,X,F,D,X) 

where  N,  X,  F,  and  D  take  the  values  returned  by  TB04AD,  and 

X  is  REAL*8,  and  is  the  point  at  which  the  value  of  the  spline  is 
required . 

I  is  INTEGER*4,  and  is  a  parameter  which  enables  the  function  to  make  use 
of  known  information  ie  if  1^0,  then  the  function  will  assume 
that  it  has  been  entered  previously  with  a  smaller  value  of  X  (this 
reduces  the  number  of  knot  intervals  to  be  searched  for  x) . 

If  I<0,  the  function  searches  the  whole  range  for  x. 


IX. 3  Function  QG02AD 


1 .  Purpose 

Given  a  cubic  spline  S(x)  defined  by  knots  x.,  function  values,  S.  =  S(x.) 
and  derivative  values  g.  =  dS(x.)/dx,  i=l,2,...,n,  n£2,  and  given 
integration  limits  a  and  bxthe  routine  will  evaluate 


b 

Q  =  /  S(x)dx 
a 


The  routine  is  not  restricted  to  splines,  and  may  be  used  with  any 
piecewise  cubic  with  continuous  first  derivative  defined  by  its  values  and 
first  derivative  values  at  the  joins. 

2.  Argument  list 

DOUBLE  PRECISION  FUNCTION  QG02AD(A,B,N,X,F,D) 

where  N,  X,  F,  and  D  take  the  values  returned  by  TB04AD,  and 

A,B  are  REAL*8,  and  must  be  set  by  the  user  to  the  lower  and  upper  limits 
of  region  of  integration  respectively.  If  either  A  or  B  is  outside 
the  range  of  the  knots,  the  integral  is  calculated  on  the  assumption 
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that  S(x)=0  for  x<x^  and  x>x  .  If  A>B  the  sign  of  the 
integral  is  appropriately  modified. 


IX. 4  Subroutine  QG05AD 
1 .  Purpose 

To  evaluate  the  definite  integral 


b 

I(a,b;  f)  =  /  f(t)dt 
a 


where  f(t)  is  defined  by  a  user-supplied  subroutine,  to  specified 
accuracy . 

2.  Argument  list 

SUBROUTINE  QA05AD ( VAL , FUNC , A , B , AERR , RERR , LEVEL , ERROR , IFLAG) 

VAL  is  a  REAL*8  variable  in  which  the  value  of  I(a,b;  f)  is  returned. 

FUNC  is  the  name  of  a  user-supplied  function  subprogramme  to  evaluate  f(t) 
for  any  value  of  t  in  the  range  [A,B).  Its  only  argument  is  the 
value  of  t,  which  is  passed  in  REAL*8;  the  value  is  also  returned  in 
REAL*8. 

A,B  are  REAL*8  input  arguments  giving  the  lower  and  upper  limits  of 
integration. 

AERR  are  REAL*8  input  arguments  specifying  the  required  accuracy  as  RERR 

max (AERR, RERR-ABS (VAL)) . 

LEVEL  is  an  INTEGER*4  input  argument  specifying  what  printed  output  (if 

any)  is  required  from  the  subroutine.  LEVEL=1  specifies  no  output. 

ERROR  is  a  REAL*8  output  argument  in  which  the  subroutine  returns  an 
estimated  upper  bound  on  the  error  achieved. 

IFLAG  is  an  INTEGER*4  output  argument  indicating  what  difficulties  were 
encountered,  if  any.  Values  of  IFLAG  up  to  3  indicate  that  VAL 
should  be  accepted  as  the  result  (subject  to  value  of  ERROR); 
the  result  may  well  be  good  even  if  IFLAG  is  4  or  5. 


lA 


1 
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APPENDIX  X  CALCOMP  PLOTTER  ROUTINES  USED  BY  THE  DATA  REDUCTION  PROGRAMME 

The  data  reduction  programme  calls  on  various  CALCOMP  plotter  subroutines.  The 
purpose  and  arguments  of  these  are  described  below.  All  dimensions  given  are 
in  inches. 

X.l  Subroutine  PLOTS 

1 .  Purpose 

Initiates  plotting  of  a  logical  segment. 

2.  Argument  list 

CALL  PLOTS ( I BUF , NO , NUM ) 

IBUF,  NO,  and  NUM  are  INTEGER*4  dummy  variables. 

X.2  Subroutine  ENDPLT 

1 .  Purpose 

Ends  a  logical  segment.  The  new  segment  will  start  with  the  pen  in  the 
standard  initial  position  and  at  the  new  logical  origin. 

2.  Argument  list 

CALL  ENDPLT 

X.!J  Subroutine  ENDGR 

1 .  Purpose 

Ends  plotting  on  the  current  device. 

2.  Argument  list 

CALL  ENDGR 

X.4  Subroutine  PLOT 

1 .  Purpose 

Moves  the  pen  to  a  specified  position. 
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2.  Argument  list 

CALL  PLOT(X,Y,NPEN) 

X,Y  are  REAL--4 ,  and  contain  the  logical  coordinates  of  the  next  pen 
position. 

NPEN  i;  INTEGER*4  and  is  the  signed  sum  of  two  positive  integers,  N  and  K. 
Il  N=2,  the  pen  is  down  during  movement. 

Il  N=3,  the  pen  is  down  during  movement. 

Normally  K=0. 

If  K=20,  the  logical  origin  is  re-located  at  (X,Y)  after  movement. 

If  NPEN=999,  the  plot  file  is  ended. 


X.5  Subroutine  SYMBOL 

1 .  Purpose 

Writes  symbols  and  characters. 

2.  Argument  list 

CALL  SYMBOL (X , Y.SIZE, TEXT , THETA , N) 

X,Y  are  REAL*4 ,  and  define  the  logical  coordinates  of  the  lower  left-hand 
corner  of  the  first  character  to  be  drawn. 

SIZE  is  REAL*4  and  specifies  the  character  height.  It  is  recommended  that 
this  be  a  multiple  of  0.07. 

TEXT  may  be  either  REAL*4  or  INTEGER*4,  and  contains  a  string  of  characters 
to  be  written. 

THETA  is  REAL*4,  and  is  the  angle  (in  degrees)  of  orientation  of  the 
character  string  in  TEXT  with  respect  to  the  positive  x  axis. 

N  is  INTEGER*4,  and  is  the  number  of  characters  to  be  written. 


X.6  Subroutine  NUMBER 

1 .  Purpose 

Writes  a  floating  point  number. 

2.  Argument  list 

CALL  NUMBER(X,Y, SIZE, FPN, THETA, N) 


X,  Y,  SIZE  and  THETA  are  as  for  SYMBOL  (see  above). 

FPN  is  REAL*4,  and  is  the  floating  point  number  to  be  written. 
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N  is  INTEGER'^,  and  specifies  the  precision  required.  FPN  is  rounded  to 
this  precision. 

If  N  is  positive,  it  is  the  number  of  digits  to  the  right  of  the 
decimal  point. 

If  N  is  zero,  the  integer  portion  and  decimal  point  are  drawn. 

If  N  is  -1,  only  the  integer  portion  is  drawn. 

If  N  is  less  than  -1,  ±N±-1  digits  are  truncated  from  the  integer 
portion . 


X.7  Subroutine  LINE 


1 .  Purpose 

Draws  a  line  througii  a  set  of  coordinates. 

2.  Argument  List 

CALL  LI NE (X , Y , N , INC , LINTYP , INTEQ) 

X,Y  are  REAL*4  arrays,  and  contain  the  x-  and  y-coordinates  of  the  data 
points  in  the  first  N  locations.  X(N+2)  contains  a  scaling  factor, 
DELTAV  (number  of  units  per  inch  on  the  x-axis),  and  X(N+1)  contains 
FIRSTV,  a  suitably  rounded  approximation  to  the  first  value  to  be 
plotted.  The  contents  of  Y(N+1)  and  Y(N+2)  are  similarly  defined. 

N  is  INTEGER*^,  and  is  the  number  of  data  points  in  X  or  Y. 

INC  is  INTEGER*4,  and  is  the  increment  of  successive  data  points  in  the 

X  and  Y  arrays  to  be  plotted. 

LINTYP  is  1NTEGKR*4,  and  controls  the  type  of  line; 

tl.INTYPi  is  the  t requency  of  plotted  symbols.  eg  if  ±LINTYP±=4 

a  symbol  (specified  by  INTEQ)  is  plotted  every  4th  data  point. 
If  LINTYP=0,  i  symbols  are  plotted. 

If  LINTYP  is  j  itive,  a  straight  line  connects  the  data  points. 

If  LINTYP  is  negative,  only  symbols  are  drawn. 

INTEQ  is  INTEGER*4,  and  specifies  a  special  symbol  to  be  drawn  at  the  data 
points . 


TEMPERATURE  (DEG.C).  temperature  IDEG.C)  . 
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Figure  2.  Temperature-salinity  diagrams  for  the  South-West  Pacific 
Ocean  constructed  from  historical  Nansen  cast  data.  The 
solid  curves  represent  a  least-squares  cubic  spline  fit  to 
the  data  (see  text) 
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Figure  4.  Profiles  of  versus  depth  constructed  using  the  temperature 

profiles  of  figure  3  and  corresponding  salinity  data.  The 
solid  curves  represent  a  least-squares  fit  of  equation  (2) 
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Figure  4(Contd.) 
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00  TWO  DEGREE 
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Figure  9.  Modified  Canadian  10°  square  subdivisions.  All  subdivisions 
of  10°  squares  are  numbered  identically  regardless  of  quadrant 
of  the  globe 
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The  expendable  bathythermograph  is  a  convenient  and  inexpensive  tool 
for  the  measurement  of  ocean  temperature  structure.  It  does  not 
measure  salinity,  but  in  many  areas  of  the  ocean  it  is  reasonable 
to  use  standard  temperature-salinity  (T-S)  curves  to  infer  salinity 
from  temperature  measurements.  T-S  curves  constructed  from 
historical  Nansen  cast  data  are  presented  for  10°  squares  in  the 
South-Western  Pacific  Ocean,  covering  the  area  from  15° S  to  45° S, 
and  from  150°  F;  to  180°  E.  With  a  knowledge  of  both  temperature  and 
salinity  structures  a  wide  range  of  other  parameters  may  be 
estimated,  the  most  important  of  these  being  density,  dynamic  height 
and  sound  speed.  Relationships  for  the  computation  of  these  and 
other  properties  of  sea  water  are  discussed,  and  computer  programmes 
presented  to  calculate  profiles  of  the  relevant  parameters  from  XBT 
data. 
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